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1..  STUDIES  RELATED  TO  THE  INVERSION  OF  THE  RADIATIVE 
TRANSFER  EQUATION:  VTPR  TEMPERATURE  SOUNDINGS  AND 
ATMOSPHERIC  COMPOSITION 


Contributors  to  the  effort  reported  in  this  section 
were  Toran  Hostbjor  (Section  1.4)  and  N.  Goravanchi 
(Section  1.3). 


This  section  reports  on  the  contractual  efforts  under 
Requirements  3.2.5,  3.2.6  and  part  of  3.2.3.  All  of  the  sections  relate 
to  probable  errors  in  VTPR  sounded  temperature  profiles  and  have  in  common 
the  inversion  of  the  radiative  transfer  equation.  The  underlying 
philosophy  of  our  effort  under  this  section  is  to  develop  techniques 
that  will  provide  accurate  VTPR  temperature  profiles  under  cloudy  conditions 
over  land  masses.  This  requires  techniques  that  will  not  only  adequately 
correct  for  c’oud  conditions  but  will  also  compensate  for  atmospheric 
composition  fluctuations. 

Under  requirement  3.2.6  we  have  undertaken  (and  completed  except 
for  checking  for  possible  errors)  a study  to  determine  the  improvement 
in  VTPR  temperature  profiles  if  SMS  infrared  type  cloud  data  were  incorporated 
into  the  retreival  scheme.  Our  preliminary  result  indicates  an  improvement 
by  at  least  a factor  of  2 (higher  altitude  temperatures  have  greater 
improvement) . 

In  either  evaluating  errors  in  VTPR  temperatures  or  in  evaluating 
possible  errors  in  composition  retrievals,  one  must  know  the  relationship 
between  a radiosonde  temperature  profile  and  an  area  average  (VTPR) 
temperature  profile.  This  evaluation  (i.e.  the  determination  of  the 
extensiveness  of  temperature  differences  for  an  error  to  exist)  is  basic 
to  the  evaluation  of  possible  errors.  This  study,  in  two  parts,  is  almost 
complete.  It  is  reported  on  in  section  1.2. 

Efforts  to  evaluate  VTPR  temperature  errors  resulting  from  probable 
atmospheric  composition  errors,  3.2.3,  have  been  directed  toward  the 
evaluation  of  possible  C0?  fluctuations.  The  approach  has  been  to  combine 


this  study  with  the  effort  under  3.2.5  to  obtain  evidence,  if  any,  of 
CO^  deviations.  This  ccald  then  be  combined  with  a previous  study  (see 
Appendix)  to  evaluate  probable  errors.  The  preliminary  result  of  effort 
3.2.5  is  tl  at  there  is  a variation  of  the  high  altitude  CO ^ content. 

The  result  cannot  be  applied  to  lower  altitudes  because  we  did  not  have 
infrared  cloud  data  for  the  times  coincident  with  the  inverted  VTPR  data. 
The  result  of  this  general  study  is  reported  in  section  1.4. 

In  summary,  the  results  of  these  studies  indicate: 

1.  Radiomctrical ly  derived  C0o  atmospheric  content  at  high  altitudes 
(25  km)  is  less  than  320  ppm  during  daylight  hours. 

2.  The  radiometric  technique  can  be  used  to  determine  atmospheric  1^0 
vapor  content,  provided  the  cloud  problem  can  be  solved. 

3.  If  SMS  infrared  data  can  be  combined  into  VTPR  retrieval  techniques 
substantial  reductions  in  temperature  errors  will  result. 

The  major  direction  of  future  work  under  these  and  related  efforts 
should  be  directed  toward: 

1.  Incorporating  (and  evaluating)  cloud  corrections  in  a watervapor 
determination. program.  This  program  should  be  based  upon  the  CO^ 
iterative  technique. 

2.  Devising. an  adequate  C07  transmittance  model  (preferably  a Smith 
Polynomial  of  at  least  2nd  or  3rd  order  in  composition)  for  use  in 
the  CO^  retrieval  effort.  The  cloud  correction  scheme  used  for  1^0 
should  be  used  in  this  effort. 

3.  Examine  the  possibility  of  incorporating  the  above  two  studies  into 
the  VTPR  retrieval  techniques. 
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1.2  STUDIES  ON  THE  EXPECTED  DIFFERENCES  BETWEEN  RADIOSONDE 

TEMPERATURE  PROFILES  AND  VTPR  SOUNDED  TEMPERATURE  PROFILES 


One  of  the  difficulties  in  evaluating  the  accuracy  of  radiometrical ly 
determined  temperature  profiles  is  the  uncertainty  in  their  relationship 
with  radiosonde  measurements  made  in  or  near  the  viewed  region.  This 
difficulty  is  further  complicated  if  the  radiosonde  measurement  is 
temporally  displaced. 

A study  was  made  of  this  problem  utilizing  radiosonde  data  collected 
at  WSMR  between  January  1969  and  June  1970.  The  results  of  the  first 
phase  of  the  study  (expected  difference  between  simultaneous , but  spatially 
separated,  radiosonde  and  radiometric  temperature  profiles)  was  reported 
to  the  American  Geophysical  Union  in  December  1975.  An  open  literature 
article  reporting  these  results  has  been  submitted  for  publication.  A 
copy  of  the  abstract  of  the  report  follows  this  section. 

The  second  phase  of  the  study  is  nearing  completion.  This  phase 
of  the  study  concerns  the  expected  difference  between  a spatially  and 
temporally  displaced  radiosonde  temperature  measurement  and  a VTPR  temperature 
sounding.  The  computer  program  has  been  written  and  tested  for  this  part 
of  the  study.  A preliminary  result  indicates  that  the  temperature  errors 
resulting  from  time  displacements  may  be  represented  by  a linear  relation- 
ship with  a larger  intercept  but  small  slope  (i.e.,  a large  error  for  small 
time  shifts  that  docs  not  increase  appreciably  with  moderate  time  increascs- 
up  to  three  hours). 

This  second  phase  should  be  completed  by  November  and  a written  report 
will  be  submitted  by  December  or  January. 
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ABSTRACT 


A large  sample  of  temperature  observations  was  statistically  analyzed 
to  estimate  horizontal  temperature  variability  as  a function  of  dis- 
tance. These  estimates  were  determined  at  several  altitudes  from 
the  surface  to  15  km  and  are  applicable  to  horizontal  distances  up 
to  175  km. 

The  results  were  used  to  assess  the  amount  of  disagreement  one 
should  expect  when  comparing  radiosonde  temperature  measurements 
with  corresponding  measurements  derived  from  satellite  radiometric 
observations.  The  conclusion  was  that  horizontal  temperature 
variations  over  the  radiometrically  observed  area  contribute  approx- 
imately IK  root  mean  square  (rms)  disagreement  between  such  compari- 
sons. When  the  separation  distance  between  the  radiosonde  and 
satellite  observation  approaches  200  km,  rms  differences  of  greater 
than  2K  are  to  be  expected. 


1.3  OPTIMIZED  STUDY  OF  USING  SMS  INFRARED  CLOUD 
DATA  IN  VTPR  TEMPERATURE  SOUNDINGS 


This  study  has  been  undertaken  to  determine  the  extent  of  possible 
improvements  in  VTPR  temperature  soundings  if  SMS  infrared  cloud  data 
were  utilized  in  the  inversion  programs.  It  is  motivated  by  the  fact 
that  the  SMS  data  presently  contains  the  cloud  data  needed  for  the 
VTPR  retrievals  and  secondarily  that  a VTPR  instrument  may  some  day  be 
included  on  the  geostationary  satellites.  Explicit  in  the  study  are  the 
assumptions  that  automatic  computer  processes  will  be  devised  for  the 
SMS  data  that  will  solve  the  registration  problem  and  provide  the 
necessary  cloud  data  for  the  VTPR  retrieval. 

It  is  assumed  that  the  basic  raw  data  available  from  the  SMS 
infrared  processing  will  be  an  average  percent  cloud  cover  and  an 
average  cloud  top  temperature  for  each  VTPR  viewed  spot . No  attempt  has 
been  made  in  this  preliminary  study  to  evaluate  the  slant  path  effect  on 
determining  cloud  percentages  and  cloud  top  temperatures  nor  has  any 
attempt  been  made  to  evaluate  differences  in  the  effective  cloud  top 
temperatures  for  the  two  instruments.  Grey  scale  registration  for  the 
SMS  data  also  has  not  been  considered.  These  problems  are  suggested  as 
possible  extensions  of  the  present  work. 

In  order  to  examine  the  possible  improvements,  a synthetic  study 
was  undertaken.  The  Atmospheric  Sciences  Laboratory  temperature  retrieval 
program,  SATFA,  was  used  to  determine  the  temperature  profile  for  each 
spot  of  a 7 X 7 VTPR  array.  Three  different  sets  of  meteorologically 
consistant  cloud  top  temperatures  and  cloud  cover  percentages  for  each 
spot  of  the  7X7  array  were  generated.  The  cloud  top  data,  the  temperature 
profiles  and  an  assumed  watervapor  profile  (not  consistant  with  the  cloud  data) 
were  used  to  compute  synthetic  radiances  for  the  VTPR. 


The  single  watervapor  profile,  the  radiances,  the  cloud  percentages 
and  the  cloud  temperatures  for  each  spot  were  used  as  input  in  a modified 
retrieval  program.  In  this  modified  program  clear  channel  radiances 
were  not  computed.  The  only  modification  made  was  in  the  computation 
of  trial  radiances  in  the  iterative  scheme;  weighting  function  were  not 
modified.  The  trial  radiances  were  computed  by  numerically  integrating 
the  transfer  equation  in  two  parts:  the  sum  for  pressure  levels  above 
the  cloud  top  was  left  unmodified;  the  sum  below  the  cloud  top  was 
modified  by  the  factor  (100  - Cloud  Percentage)/100.  The  cloud  height 
pressure  was  determined  each  cycle  using  the  latest  computed 
temperature  profile.  The  modified  subroutine  TNVFRT  of  program  SATFA 
is  shown  in  Table  I. 

Results : 

The  temperatures  retrieved  using  the  above  technique  were  compared  to 
the  actual  temperatures  used  in  the  radiance  calculations.  Typically 
the  RMS  temperature  difference  for  pressure  levels  above  300  mb  was 
less  than  1°K.  RMS  temperature  differences  for  the  entire  profile 
were  typically  in  the  range  1.6°K.  to  1.8°K. 

A comparison  was  made  of  the  retrieved  temperature  from  the  new 
scheme  with  those  determined  using  program  SATFA.  The  RMS  temperature 
errors  in  the  new  scheme  were  substantially  (approximately  1/2) 
smaller  than  the  RMS  temperature  errors  using  program  SATFA.  A typical 
computer  output  for  this  comparison  is  shown  in  Table  II. 
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The  comparison  with  program  SATFA  indicated  a possible  error  in  the 
study.  Examination  of  the  SATFA  result  indicated  errors  upward  toward 
5°K.  and  in  some  cases  even  greater.  Jlie  "SATFA"  RMS  profile  errors 
were  between  3°K.  and  4°K.  The  author  feels  that  these  errors  are  excessive 
and  error  search  has  been  undertaken.  As  of  this  writing  no  apparent 
error  has  been  uncovered;  howeverjall  data  is  being  rerun  since  it  is 
the  subject  of  the  latter  author's  thesis.  Upon  completion  of  the 
rcrun^a  special  report  documenting  the  results  will  be  prepared. 

Even  though  an  error  is  suspected,  we  do  not  believe  that  the  general 
character  of  the  result  is  wrong.  The  revised  technique  produces  excep- 
tionally small  temperature  error  for  altitudes  greater  than  300  mb. 

This  result  should  not  change.  While  lower  altitude  retric/ed  temperatures 
contain  less  error,  the  improvement  is  less  dramatic. 
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1.4.1  SUMMARY 

NOAA-2  satellite  radiance  measurements  and  correlated  radiosonde 
temperature  and  dew-point  temperature  profiles  for  the  WSMR  Region  were  used 
to  determine  the  variability  of  atmospheric  CO.,,  for  different  altitudes, 
by  asymptotic  solution  of  the  Ladenberg-Reiche  equation.  Thirty-four  sets 
of  data  were  available  for  use.  Transmittances  used  in  the  retrieval  were 
corrected  using  the  measured  temperature  and  dew-point  temperature  profiles. 
Preliminary  results  indicate  a lowered  daylight  period  atmospheric  CO,, 
mixing  ratio  at  high  altitude  (as  compared  to  the  accepted  standard  mixing 
ratio).  This  appears  to  be  consistant  with  a previous  result  indicating  a 
diurnal  CO mixing  ratio  variation  (Bruce  et . al . 1975). 

An  iteration  technique,  due  to  Chahine  (1972),  was  attempted  on  the 
data  with  very  unsatisfactory  results.  Through  the  use  of  synthetic 
studies,  the  problems  in  determining  CO.,  composition  appear  to  be  the 
result  of  cloud  contamination,  inadequacy  in  the  transmittance  model  and 
the  temperature  dependence  of  the  weighing  function. 

The  study  indicates  that  the  technique  should  be  usable  for  H.,0 
determination  using  VTPR  channels  7 and  8 (and  possibly  channel  6)  if 
the  cloud  problem  can  be  satisfactorily  solved. 
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1.4.2  PERTURBATION  TECHNIQUE  RESULTS 


Basic  Formulation 

The  method  used  to  solve  the  radiative  transfer  equation  (RET)  by 
perturbation  technique  for  the  C0?  composition  profile,  q,  is  basically 
that  proposed  by  M.  T.  Chahinc  (1972).  The  outgoing  radiance  at 
frequency  v,  from  a non- scattering  atmosphere  in  local  thermodynamic 
equilibrium,  observed  from  pressure  p can  be  expressed  as 

rln  P aMT(p)] 

IV(P)  = Bv[T(p)]  - / t [q.(p),T(p)]  -d  v ■ d In  p , (1) 

An  po 

in  which  B(T)  is  the  Planck  function,  pQ  is  ground  surface  pressure  and 
t is  the  transmittance  from  pressure  level  p to  p.  It  is  noted  that  the 
transmittance  i is  a function  principally  of  the  composition  profiles, 
q^p) , of  the  several  absorbing  species  as  well  as  of  temperature,  T. 

There  is  one  such  equation  for  each  radiometer  channel.  The  usual  practice 
is  to  express  the  total  transmittance  as  the  product  of  separable 
CC^.  I^O  and  0^  transmittances , 

3 

T(q1,q2,q3,T)  = II  T (q  »T)  . (2) 

i=l 

In  this  analysis,  rather  than  solving  for  the  CO^  composition  profile 
qi(p)  we  solve  for  the  relative  CO^  profile,  a(p) , defined  by 

q(p)  - a(P)  qstd(p)  , (3) 

where  q ,(p)  is  the  CO-  composition  profile  used  in  the  standard  transmittances. 

S ill  i! 


In  order  to  solve  (1)  for  the  parameter  ci(p),  we  note  that  the  major 
contribution  to  the  outgoing  radiance  comes  from  the  atmospheric  region 
centered  about  the  peak  of  the  weighting  function;  that  is,  the  pressure 
pm  for  which  3x/3  In  p reaches  its  maximum  value.  For  a particular 
frequency  v,  solution  of  (1)  yields  a value  of  «v(Pm)  representative  of  the 
relative  composition  profile  at  pressure  pm>  Since  the  optical  depth  in  this 
region  is  on  the  order  of  1 or  greater,  strong  line  absorption  is  valid 
(Kaplan,  1953) . The  parameter  can  be  related  to  t standard  through 
asymptotic  solution  of  the  Ladenberg-Reiche  equation  (Chahine,  1972)  , 


x (ci  q ,,T)  ~ T (q  ..J) 
C.O2  v std  CO2  std 


In  accordance  with  Chahine  (1972)  a direct  first  order  perturbation  solution 

for  a leads  to 
v 

a - 1 - (I  - I ) 1^  \ , , in  , fifO-  6 In  pi  " . (5) 

v v v7  //.  total  co  3 In  p r{ 

In  pfl  2 ) 


in  which  is  defined  by  (2) . The  relative  CO  2 composition  in  the 

region  of  the  weighting  function  centroid  for  each  radiometer  channel  is 
obtained  from  (5)  if  one  knows  the  temperature  corrected  standard  CO, 
transmittances , 1^0  and  0^  transmittanccs  and  the  temperature  profile.  This 
approximate  solution  for  a becomes  questionable  when  it  deviates  substantially 

yf 

from  the  value  of  1,  since  the  first  neglected  term  in  the  expansion  of  x 
becomes  non-negl igib le  for  a < 0.4  or  a > 1.6.  Further,  for  a < 0.6  one 


questions  the  Ladenberg-Reiche  asymptotic  solution. 


Data 

Thirty-four  sets  of  data  correlating  temperature  and  dew-point 
radiosonde  measurements  over  WSMR,  with  NOAA-2  VTPR  soundings  were 
available  for  this  study.  All  VTPR  soundings  over  WSMR,  or  immediately 
adjacent  to  the  range,  are  being  processed  for  a CO2  perturbation 
determination.  As  of  the  close  of  the  contract  year  all  of  the  available 
data  had  not  been  processed.  A reasonable  amount  of  time  is  required  to 
select  the  specific  VTPR  soundings  to  be  processed  from  the  49  blocks  in 
each  set.  Each  sounding  location  must  be  plotted  on  a map  of  the 
WSMR  region. 

Results 

Figure  IV- 1 through  1V-7  illustrates  the  typical  results  obtained 
using  this  technique.  Figure  IV-1  and  IV-2  show  the  radiosonde  water 
vapor  and  temperature  profiles  for  the  particular  sounding  set.  The  next 
five  figures  show  the  inferred  CO^  mixing  ratios  for  each  of  the  five 
VTPR  soundings  on  or  adjacent  to  WSMR.  The  vertical  axis  for  each  plot 
is  the  standard  VTPR  transmittance-pressure  level,  numbering  1 at  the  top 
of  the  atmosphere  and  100  at  the  bottom.  (The  number  "101"  is  an  artifact  of 
the  plotter.)  Each  of  the  C07  profile  plots  show  C0^  for  the  measured 
temperature  and  watervapor  • circles ; and  the  probable  limits  of  accuracy 
assuming  10%  watervapor  fluctuation  and  +_  1°K.  in  the  temperature  measurement 
Discussion 

Comparison  of  the  results  with  the  temperature  profile  indicates  that  th 
CO2  determination  for  pressure  level  55  (VTPR  channel  3)  is  invalid--the 
temperature  is  essentially  isothermal  in  this  region.  The  determination 
for  pressure  level  45  (VTPR  channel  2)  is  very  suspect  for  similar  reasons. 


1.4.3  ITERATION  TECHNIQUE  RESULTS 


Basic  Formulation 

As  originally  formulated  by  Smith  (1970)  and  extended  by  Chahine  (1972) 
equation  1 may  be  solved  for  the  scaling  constant  a,  equation  3,  using 
an  iterative  technique.  Since  a separate  ok  is  determined  for  each  VTPR 
channel,  the  technique  requires  convergence  on  an  atmospheric  profile 


defined  by 


(6) 


where 


cu  = scaling  constant  for  the  j ^ VTPR  channel 

■h  Vi 

P.  = i VTPR  pressure  level 

t h 

W (Vj , P^)  = weighting  function  for  the  j VTPR  channel  at  pressure  level  P . 

Determination  of  the  weighting  function  requires  a model  for  the  CO^ 
transmittance:  the  n^1  iteration  weighting  function  is  given  by 
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in  which  B(vj,  P^)  is  the  Planck  function  for  the  j1'1  VTPR  channel  frequency 
at  pressure  level  P. . It  is  through  Planck  term  that  the  temperature  profile 
enters  the  problem.  It  is  obvious  that  for  an  isothermal  temperature  region, 
no  composition  information  can  be  obtained. 

For  the  iteration  technique  the  NOAA-2  zenith  angle  correction  model 
was  utilized.  This  model  is  inadequate  since  only  the  linear  terms  in  a 
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of  a Smith  Polynomial  model  are  Used: 
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in  which  for  the  jt*1  channel. 
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Using  this  model  the  first  term  in  the  weighting  function  for  the  j 
VTpR  channel  and  pressure  level  P^  is  given  by 

: rt  P.)  Oo) 
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An  iterative  computer  scheme  based  upon  the  Chahine  technique,  which 
is  similar  in  nature  to  the  temperature  retrieval  iterative  schemes,  was 
written  for  the  University  IBM  computer. 

Computer  Results  and  Difficulties 

Excessive  convergence  difficulties  were  encountered  when  the  WSMR  data 
was  utilized.  The  primary  suspected  cause  for  this  problem  (outside  of 
programming  error)  was  the  transmittance  model . To  test  the  adequacy  of 
the  program  and  the  transmittance  model,  synthetic  data  was  prepared: 
transmittances  for  various  CO,,  profiles  were  computed  using  the  model  and 
in  turn  used  to  compute  synthetic  radiances  for  the  observed  temperature 
and  watervapor  for  one  particular  set  of  data.  These  synthetic  radiances, 
the  temperature  profile  and  watervapor  profile  were  used  as  input  for  the 
iteration  program.  While  convergence  was  slow^it  was  obtained- and  the 
resulting  computed  CO^  profiles  for  all  but  one  case  were  within  tolerable 
limits  of  the  correct  value.  Computer  output  for  the  convergent  cases 
arc  shown  in  Figures  IV-8  through  IV-17.  The  assumed  temperature  profiles 
and  watervapor  profiles  used  for  the  calculation  are  shown  in  Figures 
IV-8  and  IV-9.  The  remaining  figures  in  the  set  are  ordered  so  that  the 
CO2  profile  used  in  the  radiance  calculation  is  followed  by  the  retrieved 
CO^  profile. 

Figure  I V- 1 8 shows  the  final  weighting  functions  for  case  4.  Comparison 
of  the  retreived  profile  with  the  actual  profile  for  case  4 (Figures  IV-17 
and  IV- 16)  shows  a discrepancy  at  both  the  top  and  the  bottom  of  the 
atmosphere.  This  is  explained  by  the  fact  that  the  weighting  functions 
fall  off  to  zero  at  their  pressure  levels.  Figure  1 V - 1 8 also  indicates  the 


p 

effect  of  the  isothermal  region  of  the  temperature  profile:  little  or 
no  information  can  be  obtained  about  the  region  near  the  tropopause.  We 
believe  that  this  effect  may  be  the  cause  of  the  non-convergence  in  case  5 
(Figures  IV-19  and  IV-20) : too  little  CO^  is  determined  in  the  tropopause 
region  and  too  much  in  the  lower  atmosphere.  This  case  also  probably 
indicates  a difficulty  with  the  model;  second  order  terms  in  a are  needed. 

Discussion 

The  synthetic  study  indicates  that  the  computer  algorithms  are 
programmed  correctly.  We  feel  that  the  lack  of  convergence  using  actual 
data  can  be  traced  to  either  a)  inadequate  transmittance  model  or  b)  cloud 
contamination,  or  a combination  of  the  two.  The  former  of  the  two  is 
clearly  indicated  by  examination  of  figure  IV-18:  essentially  only  two 
regions  of  the  atmosphere  can  be  sensed  for  CO,,  using  the  inadequate 
transmittance  model. 

At  present  we  are  working  on  obtaining  a better  C0o  model.  With  this 
anticipated  improvement  and  assuming  that  the  clear  channel  radiances 
can  be  determined,  we  believe  that  it  is  possible  to  obtain  good  CO^ 
profiles.  Lvcn  without  cloud  corrections,  it  appears  that  upper  channel 
CO2  values  can  be  obtained. 


i 
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1.4.4  RECOMMENDATIONS  FOR  FURTHER  WORK 


So  that  this  effort  can  be  brought  to  a successful  conclusion,  we 
recommend  that  the  following  be  undertaken. 

a)  A more  adequate  CO^  transmittance  model  for  NOAA-2  transmittances 
should  be  constructed.  We  believe  that  a Smith  type  polynomial  containing 
at  least  second  order  terms  in  composition  (and  preferably  third  order) 

is  needed.  Work  on  this  will  continue. 

b)  Cloud  correction  techniques  should  be  applied  to  the  radiances 
used  for  the  CO  2 retrievals.  At  this  time  it  is  not  clear  which  of  the 
several  cloud  correction  techniques  reported  elsewhere  herein  will  be 
the  most  effective. 

c)  The  iteration  technique  should  be  utilized  to  determine  H^O 
content  in  the  temperature  retrieval  program.  A preliminary  effort 

is  underway  to  demonstrate  that  channels  7 and  8 or  channels  6,  7 and  8 
can  be  used  for  this  purpose.  Of  great  interest  will  be  the  effect  of 
water  normalization  using  a trial  temperature  on  the  retrieved  temperature 
profile . 

d)  Cloud  correction  techniques  should  be  incorporated  into  the  1^0 


retrieval  subroutine. 
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Figure  IV-4 
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Figure  IV-6 
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Figure  IV-9.  Water  vapor  profile  use  for  synthetic  study. 
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Figure  1V-14.  Actual  C0_  profile  Case  3. 
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Figure  IV- 19.  Actual  CO-  profile  for  Case  5 


Section  3.2.3  of  the  contract  requirements  called  for  the  production 
of  homogeneous  path  transmittance  calculations  necessary  for  the  transmittance 
model  work.  During  the  course  of  the  contract  year,  400  different,  filter 
convoluted,  homogeneous  path  Carbon  Dioxide  transmittances  were  computed 
for  each  of  the  six  VTPR,  NOAA-2  CO.,  filter  channels.  These  transmittances 
were  computed  for  different  CO.,  gas  amounts  for  standard  temperature  and 
pressure.  In  addition  100  such  ozone  and  100  line  contribution  water 
vapor  transmittances  were  computed  for  each  of  the  first  six  NOAA-2  filter 
channels . 

Section  3.2.4  called  for  an  examination  of  errors  resulting  from  the 
standard  technique  of  using  separate  C0?,  0^  and  H^O  transmittances  in 
the  VTPR  retreival  program.  A preliminary  line-by-line  filter  convoluted, 
transmittance  calculation,  for  an  atmosphere  of  10  layers,  was  made 
utilizing  H^O  and  CO^  lines  for  VTPR  channels  5 and  6.  Ihese  two 
transmittances  for  the  H^O-CO^  model  atmosphere  were  compared  to  the 
product  of  separate  filter  convoluted  llo0  and  CO.,  transmittances  for  the 
same  10  layer  atmosphere.  For  channel  5 the  combined  transmittance  was 
greater  than  the  product  while  for  channel  (>,  the  exact  opposite  result 
was  observed.  As  far  as  this  writer  is  concerned,  there  appears  to  be  no 
physical  or  mathematical  reason  (other  than  error)  to  explain  this  result. 
After  a considerable  error  search  in  the  computer  programs  it  was  determined 
that  the  random  round  off  error  in  the  line-by-line  computational  technique 
was  the  only  logical  source  of  error.  A test  was  made  and  the  cumulative 
round  off  error  over  a 50  wave  number  band  was  determined  to  vary  upward 
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to  .2  wave  numbers.  A new  program  utilizing  a different  computational 


scheme  for  slant  path  calculations  with  considerably  increased  double 


precision  summations  has  been  written  and  tested.  The  output  of  this 


program  is  a magnetic  tape  of  high  resolution  layered  transmittance  data 


The  second  phase  program  of  the  new  procedure  is  in  the  process  of  being 


This  second  program  reads  the  stored  data,  modifies  for  composition 


variation,  computes  the  product  transmittances  through  the  atmosphere  and 


convolutes  the  data  with  the  VTPR  filters.  When  these  programs  are  used 


the  entire  CO.  band  can  be  computed  with  an  assured  accuracy  that  each 


point  is  accurate  in  wavenumber  to  .01  cm  (It  will  require  more  than 


one  magnetic  tape  to  compute  an  entire  atmosphere  with  higher  resolution.) 


When  the  second  program  is  completed,  the  desired  calculations  can 
be  made  with  assurance  that  the  random  displacement  of  the  computed  line 


data  will  be  eliminated 


2.2.1 


SUMMARY  OP  RADIATIVE  TRANSFER  EQUATION  PROBLEM 


The  fundamental  assumptions  basic  to  the  standard  techniques  used 


in  VTPR  inversions  is  that  the  atmosphere  can  be  represented  as  plane 


parallel,  in  local  thermodynamic  equilibrium  and  free  of  scattering  agents, 


i.e.  no  multiple  scattering.  For  such  an  atmosphere  the  radiative  transfer 


equation  (RTE)  has  the  solution  (for  monochromatic  radiation  at  frequency  v) 

i 

l(v>,0)  = * (BLv/T(p)]  , 1 

in  which  the  subscript  o indicates  the  lower  boundary,  and 


p = pressure. 


B[v,T(p)]  = Planck  Function, 
0 = Zenith  Angle, 


(P  ,j 

- exp  J-  I k(v,  p‘)  ^(p’J  Sec  O ‘ip 


In  equation  2,  k(v,  p ) is  the  absorption  coef  at  frequency  v for  gas  q(p  ). 


The  difficulty  in  VTPR  retrieval  is  that  the  radiometer  responds  over 


a broad  band  of  frequencies.  Rigorously  equation  1 should  be  convoluted 


with  the  filter  band  pass  of  the  VTPR,  i.e.  integrated  over  frequency. 


If  this  were  done  the  problem  is  unsolvable.  Therefore,  the  approximation 


is  made  that  the  plank  function  can  be  represented  by  an  average  value: 


its  response  at  a single  wavenumber,  \k  , in  the  bandpass  interval.  With 


this  approximation  the  pressure  and  frequency  integrations  are  interchanged 


and  weighted  transmittances  are  computed  to  replace  the  monochromatic 


transmittances : 
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in  which  fj(v)  is  the  bandpass  function  (filter  transmittance)  and 


and  v 2 are  the  points  beyond  which  f j (v)  = 0.  The  format  of  equation  1 


is  unchanged  except  all  terms  now  refer  to  frequency  averaged  values. 

It  is  apparent  that  some  error  has  been  introduced  into  the  solution 
of  the  RTE  problem  even  before  mathematical  algorithms  are  applied  to  the 


solution.  The  differential  elements  dr(v.,  p,  G)  are  defined  by 
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In  order  to  solve  equation  1 in  its  modified  form,  the  transmittances 
defined  by  3 and/or  4 must  be  computed  numerically.  The  next  section 
reviews  this  problem. 
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2.2.2  GENERAL  ASPECTS  OF  TRANSMITTANCE  COMPUTATION  TECHNIQUES 


In  order  to  numerically  compute  the  incremental  monochromatic 
transmittance  at  frequency  v for  a multicomponent  gas  (composition  components 
Aq j ) at  temperature  T and  pressure  p one  must  evaluate 


A^,  . exp 


In  equation  5 all  absorption  lines,  centered  at  v^,  of  each  species  that 

contribute  to  the  absorption  at  frequency  v should  be  included  in  the 

th 

calculation.  The  absorption  coefficient,  k..,  for  the  i line  of  the 
j species  is  a function  of  the  pressure,  temperature,  frequency 
displacement  from  the  line  center  and  of  the  species  q. 

To  compute  the  filter  convoluted  transmittance  for  a particular 
homogeneous  path,  equation  S must  be  summed  over  all  frequency  points, 


A A , within  the  filter  bandpass: 
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in  which  the  frequency  points  are  equally  spaced. 

For  the  computation  of  nonhomogeneous  paths (atmospheric  slant  paths) , 

the  transmittance  to  each  pressure  level  at  frequency  v is  the  product 

of  all  increment  transmittances  for  the  various  layers: 

A 


rt)  n 
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The  convoluted  filter  transmittances  (for  the  j 
level  is  given  by 


assuming  equally  spaced  frequency  points. 


filter)  for  each  pressure 


Am* 


The  essential  differences  between  the  computational  aspects  of  line-by-line 
programs  centers  the  line  shape  used  and  on  how  nonhomogeneous  layers  arc 
converted  to  homogeneous  layers. 

There  is  little  question  that  dopper  broadening  should  be  considered 
at  high  altitudes  and  that  the  lorentz  line  shape  is  reasonable  for 
low  altitudes.  In  the  intermediate  region,  the  voight  line  shape  is  the 
best  choice  of  the  three  on  purely  physical  grounds,  however,  this  requires 
substantial  increases  in  computer  time.  In  addition  one  may  incorporate 
collision  narrowing  effects. 

If  the  lorentz  shape  is  used,  composition  can  be  made  a function  of 
pressure  and  a closed  analytical  form  can  be  obtained  for  the  integral 
over  pressure  (assuming  constant  temperature) . Other  choices  to  convert 
nonhomogeneous  layers  to  homogeneous  layers  is  to  use  median  values  of 
P and  T or  convert  to  a pressure  weighted  mean  temperature  for  each  layer. 
There  are  several  non>-computational  technique  options  available.  These 
include  the  selection  of  lines  to  be  included  in  the  calculation,  whose 
set  of  line  parameters  will  be  used,  the  number  of  layers  in  a vertical 
path,  the  extensiveness  of  wing  contributions  and  the  frequency  density 
of  the  calculations.  The  availability  of  data  usually  determines  whose 
line  parameters  are  used.  ITte  other  choices  arc  many  times  determined  by 
computer  time  and  space  considerations.  For  instance,  the  choice  of  low 
excitation  energy  levels  for  low  temperature  altitudes  and  high  line 
strength  for  hotter  altitudes  is  usually  compromised  by  the  availability 
of  program  storage  space. 


2.2.3 


DATA  AVAILABLE  FOR  CALCULATIONS  AND  ERROR  SOURCES 


> 


i 


All  calculations  made  use  of  the  "McClatchey  Data  Tape"  for  line 
parameter  information.  It  is  apparent  that  some  minor  errors  and  omissions 
exist  in  the  data  because  even  the  most  accurate  calculations  for  homo- 
geneous paths  show  some  differences  in  the  comparison  to  equivalent 
experimental  data.  Since  this  data  is  periodically  upgraded,  the  errors 
and/or  omissions  will  reduce  with  time. 

Line-by-line  calculations  for  the  same  atmospheric  path  computed  by 
different  agencies  can  vary  by  up  to  5%  (as  recently  reported  at  the 
Workshop  on  Satellite  Atmospheric  Soundings) . This  does  not  seem 
unreasonable  when  one  considers  the  possible  choices  inherent  in  the 
calculations.  It  does,  however,  indicate  considerable  need  for  a more 
careful  analysis  of  possible  error  sources  inherent  in  the  computational 
techniques.  The  main  difficulty  in  evaluating  the  errors  is  the  lack  of 
usable  slant  path  experimental  data;  observations  over  slant  paths  in  which 
one  has  good  knowledge  of  gas  composition,  temperature  and  pressure. 

Apparently  the  only  slant  path  data  availableare  the  VTPR  radiance 
measurements  obtained  with  simultaneous  radiosonde  measurements.  An 
investigation  of  transmittance  errors  utilizing  this  data  will  by 
necessity,  be  a statistical  study  and  also  will  be  time  consuming.  For 
the  present  we  can  see  no  better  or  practical  approach  to  the  error  problem. 
The  new  slant  path  computation  scheme  we  have  devised  will  make  a study 
of  this  type  somewhat  easier  since  composition  modifications  can  be  handled 
much  easier. 
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2.2.4 


PHILOSOPHY  OF  CALCULATION  APPROACH 


Until  a detailed  study  of  the  error  sources  is  completed  we  believe 
that  working  transmittances  should  be  based  upon  line-by-line  calculations 
using  only  doppler  and  lorentz  profiles.  The  other  error  sources  inherent 
in  the  line-by-line  calculations  (line  parameters,  wing  contributions,  etc.) 
make  the  added  cost  of  voight  profile  computations  difficult  to  justify. 

The  present  programs  are  easily  changed  to  include  other  line  shapes. 
Modifications,  for  test  purposes,  will  be  made  to  include  voight  profiles 
as  specified  in  the  error  study. 


Three  versions  of  the  program  outlined  in  Appendix  1 were  used  to 
compute  CO^,  l^O  and  0^  homogeneous  transmittances  for  the  transmittance 
model  work. 

The  homogeneous  path  program  is  a modification  of  the  McClatchey 
program.  It  uses  the  lorenty  line  shape  only.  The  McClatchey  data 
are  read  from  a transfer  tape  which  contains  only  the  15u  band  data. 

The  400  CO^,  100  H^O  and  100  0^  homogeneous  path  transmittances 
for  each  filter  were  computed  for  the  following  specifications: 

1.  Lorenty  line  shape  only. 

2.  Standard  Temperature  and  Pressure-with  its  implication  on 
line  selection. 

3.  Temperature  and  pressure  corrections  for  partition  function, 
line  strength  and  line  half-width 

4.  Wing  contributions  out  to  15  cm  * 

5.  Frequency  density  of  computed  points  at  .001  cm  * intervals. 

The  program  output  not  only  gave  the  curve  of  growth  for  the  particular 
filter  but  it  also  gave  low  resolution  transmittance  data  covering  the  filter 
frequency  band.  The  low  resolution  data  was  not  needed. 

The  use  of  standard  pressure  and  temperature  in  these  calculations 
prevent  the  data  being  scaled  directly  to  low  temperature  because  the  line 
selection  would  have  been  different  for  the  low  temperature  region. 

Please  note  that  the  program  in  Appendix  2 is  still  considered  as  a 
working  program,  unfinished,  and  does  contain  a number  of  nonused  variables. 


A new  computational  scheme  has  been  devised  for  the  slant  path 
transmittance  computations.  The  program  in  Appendix  B has  been  tested 
and  will  compute  line-by-line  transmittances  for  a layered  atmosphere. 

Transmittance  data  is  generated  at  frequency  intervals  of  up  to 
.001  cm  *,  averaged  and  placed  on  magnetic  tape  at  .01  cm"'1  intervals. 

'Ihe  calculations  are  made,  for  the  most  part,  in  double  precision.  The 
present  scheme  lias  eliminated  the  random  round  off  errors  noted  in  the 
McClatchey  program  and  our  other  modifications  of  that  program. 

In  its  present  form  the  l.orentz  line  shape  only  is  utilized.  The 
equivalent  homogeneous,  temperatures  and  pressures  for  atmospheric  layers 
are  simple  median  values  for  the  layer.  Temperature  and/or  pressure  corrected 
half  widths,  line  strengths  and  partition  functions  are  included  in  the 
scheme.  The  program  will  cycle  through  a layered  atmosphere;  the  number 
of  layers  controlled  by  input  data. 

The  computational  approach  is  to  determine  transmittances  in  two 
parts.  The  program  (in  Appendix  3]  outputs  the  basic  line-by-line 
data  for  the  entire  CC^  band  region  for  each  of  the  atmospheric  layers. 

Our  preliminary  calculation  indicates  one  standard  magnetic  tape  will 
hold  a 35  layered  atmosphere  with  each  layer  containing  data  at  .01  cm  1 
intervals.  The  second  program  is  still  under  test.  It  reads  the  line-by- 
line  data  tape  (up  to  three  tapes-  3 species),  computes  the  product 
transmittances  through  the  various  layers  - one  at  a time,  and  convolutes 
the  results  with  the  various  filters. 

Iliere  are  two  main  reasons  for  this  approach: 

1.  The  line-by-line  part  of  the  computations  are  the  most  time 
consuming  part  of  the  process  - under  the  present  plan  the  data  need  be 


computed  only  once  for  each  temperature  profile. 

Composition  variations  and  different  filter  functions  can  easily 
be  incorporated  into  the  scheme  without  having  to  repeat  the  line-by-line. 

2.  By  separating  the  filter  convolution  aspect  from  the  program, 
more  computer  space  is  available  for  line  data  and  other  possible  line 
shape  routines. 

As  soon  as  the  second  program  has  been  checked  out,  a completely 
documented  report  on  the  technique  will  be  submitted. 
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1.  INTRODUCTION 

An  analysis  has  been  made  of  the  Satellite 
Temperature  Sounding  Technique  (Wark  and  Fleming, 
1966;  Smith,  1970;  Smith  et  al,  1972;  Wark  and 
Hilleary,  1969;  Wark,  1970)  to  determine  the 
extent  of  temperature  errors  arising  from  pos- 
sible atmospheric  C02  variations.  These  possi- 
ble errors  are  compared  to  probable  ozone-caused 
errors  and  to  errors  resulting  from  the  radiom- 
eter uncertainties. 

The  need  for  the  present  analysis  is  more 
than  academic.  A recent  preliminary  analysis 
of  SIRS-A  data  by  the  authors  indicates  non- 
negligible  C02  variability.  These  findings  a- 
gree  with  the  data  obtained  by  Gates  et  al  (1958), 
in  which  there  was  evidence  of  a very  nonuniform 
atmospheric  CO2  mixing  ratio.  While  these  re- 
sults are  not  experimentally  as  satisfying  as 
the  high  altitude  measurements  by  Hagemann  et  al 
(1959),  they  are  in  general  agreement  with  other 
published  evidence  of  atmospheric  CO2  variability 
in  the  lower  atmosphere.  The  measurements  re- 
ported by  Walter  Bischof,  Bert  Bolin  and  others 
(see,  for  Instance,  Bischof,  1971;  Bolin  and 
Bischof,  1970)  indicate  variations  in  time- 
averaged  C02  on  the  order  of  a few  percent. 
Woodwell  et  al  (1973)  report  substant ially 
higher  variability  near  ground  level  (125  meters). 
It  is  apparent  that  there  is  variability  in  the 
atmospheric  C02  content;  the  only  question  is 
the  amplitude  of  this  variation.  Based  upon 
these  considerations,  an  analysis  of  the  effect 
of  C02  variations  on  satellite-retrieved  tem- 
peratures seems  Justified. 

Two  methods  for  evaluating  the  temperature 
errors  were  used.  The  first  involves  develop- 
ment of  a closed  set  of  equations  for  the  ap- 
proximate temperature  error  as  a function  of 
C02  mixing  ratio  deviation  at  particular  alti- 
tudes— one  different  altitude  for  each  radiom- 
eter channel.  This  method  is  based  upon 
Chahine's  (1970)  analysis  of  the  relevant 
radiative  transfer  problem.  In  the  second 
method  the  temperature  differences  were  evalua- 
ted from  temperature  profiles  generated  by  a 
retrieval  program  in  which  the  CO  transmi t tances 


were  modified  to  reflect  varying  C02  mixing 
ratios  at  different  altitudes.  The  results  of 
the  first  method  are  in  reasonable  agreement 
with  the  results  of  the  more  accurate  second 
method  for  similar  circumstances. 

2.  CLOSED  SOLUTION 

A closed  set  of  equations  is  developed  by 
assuming  that  the  radiation  measured  by  the 
satellite  radiometer  channel  is  related  to  the 
atmospheric  parameters  at  and  above  the  centroid 
of  the  transmittance  weighting  function.  The 
differences  between  two  temperature  profiles  are 
equated  to  the  differences  in  the  two  C02  mixing 
ratios  forming  the  transmi ttances . The  atmos- 
phere is  divided  into  non-overlapping  regions, 
with  one  region  for  each  radiometer  channel 
centered  about  the  centroid  of  the  channel's 
weighting  function  (pressure  at  which  the 
weighting  function  is  maximum).  Different  C02 
mixing  ratios  are  considered  for  each  such 
region  of  the  atmosphere. 

For  simplicity  of  notation  the  channels  are 
ordered  according  to  the  altitude  of  the  maximum 
of  their  weighting  functions;  i.e.,  the  i * 1 
channel  is  the  highest.  The  atmospheric  regions 
are  Indexed  according  to  the  channel  index. 

The  measured  radiance,  I in  the  i*^  chan- 
nel, frequency  v^,  is  given  by 

*nrs  9t, 

I,  - B [T(P  )lTt(F  )+  / B,[T(P)) dlnP,  (1) 

° ° lnPQ  3enP 

where  B. (T)  is  the  Planck  function,  P is  pressure, 
and  t^(P)  is  the  transmittance.  The  subscripts 
o and  s represent  surface  and  satellite  condi- 
tions, respectively.  T(P)  is  the  temperature 
obtained  using  the  standard  C02  atmospheric  mix- 
ing ratio  (320  ppm). 

The  radiative  transfer  equation  resulting 
from  a different  C02  mixing  ratio  profile,  q (P), 
with  transmittance  t*  and  temperature  profile 
T*(P),  is 


*nPg  3t* 

= B [T*(P0)Jx*  + / lT*<P)  1 ditnP.  (2) 

lnPQ  3HnP 

The  temperature  errors,  AT,,  are  obtained  in 
terms  of  the  relative  deviation  from  standard 
in  the  CO2  mixing  ratio  by  equating  Eqs.  (1)  and 
(2)  and  expressing  the  quantities  in  the  latter 
in  terms  of  those  in  the  former. 

In  the  ic^  interval  the  mixing  ratio,  q*,  is 
related  to  the  standard  ratio,  q,t  by 

q*  * (1  + ai>  qf  (3) 

The  modified  transmi ttances , At*,  in  each 
region  are  expressed  in  terms  of  the  standard 
transmittances,  At,,  and  the  respective  u,  by 


t*  - Ax,  a*-!*. 


As  a lower  limit  approximation  for  the  effect  of 
a,,  one  could  replace  the  quantity  (l+-a,)  in 
Eq.  (4)  with  ,l+o,.  Sample  calculations  indicate 
that  the  relationship  in  Eq.  (4)  yields  more 
accurate  results  over  most  of  the  15  pm  band 
atmospheric  transmission  profiles. 

The  computations  of  the  transmittances  t*(P) 
and  the  weighting  functions  3T*/3S.nP  for  the 
mixing  ratio  q*(P)  are  shown  in  the  Appendix. 

The  Planck  functions  are  related  through  a 
first  order  expansion: 

3B, 

B,(T*(P)]  - B,[T<P)]  + AT, (P) . (5) 

3T 

Substituting  Eq.  (5)  into  (2)  results  in 


- { B, [T(P0) 1 + 

3B,[T] 

AT(P)hJ(P0)  + 

3T 

\n?B 

3B.lT] 

3tJ 

/ {B,[T(P')]+ 

AT(P'))  dUnP. 

inPQ 

3T 

SfcnP  (6) 

The  Integrals  in  F.qs.  (1)  and  (6)  are  evaluated 
according  to  Chahine's  mean  value  approach: 


Jb,(T) dlnP  = B,lT(P,)l  

31nP  3i.nP 


A1«.nP.  (7) 


The  Integral  of  Eq.  (6)  is  evaluated  at  the 
pressure  at  which  the  integrand  reaches  its 
maximum  value,  P^,  which  is  displaced  from  the 
centroid  of  the  weighting  function  in  Eq.  (1). 
The  equivalent  width  of  the  Integrands  in  de- 
fined by 

inPs  3t,  3t,  I 

A,tnP  ” [/  B,[T(P)] dlnP]/B  [T(P)] (8) 

inP„  3«nP  3tnP 


A,  = B, [T(Pj ) ] A,tnP1, 

3£nP  1 


F = B,[T(P6)]  A,tnP6, 

SdnPfi 

C,  - B,[T(P0)]  [t*(P0)  - T (P„) ] , 

3B, (T) ] „ 

D,  = t (P0)  AT0,  and 


The  subscripts  1 and  6 refer  to  the  respective 
maxima  of  the  integrands  of  Eqs.  (1)  and  (6). 

Eq.  (9)  (through  the  3?nT*/9£nP  dependence) 
shows  that  the  temperature  error  at  a particular 
level  is  related  not  only  to  the  modified  C(>2 
content  at  the  level,  but  also  to  the  CO2  varia- 
tions at  higher  levels.  In  practice  the  ground 
temperature  error  terra  must  be  determined  first, 
since  it  affects  the  error  at  other  levels 
through  the  factor  D^.  The  error  determination 
for  higher  altitudes  is  essentially  independent 
of  the  surface  term. 

For  the  purpose  of  illustrating  Eq.  (9),  tem- 
perature errors  resulting  from  uniform  CO2  atmos- 
pheric variations  were  computed  using  the  SIRS-B 
transmittances  and  the  temperature  profile  shown 
in  Figure  1.  The  resultant  temperature  errors 
for  the  C09  relative  deviations  from  standard  of 
±5%  and  ±lD%  are  shown  in  Figure  2.  In  each  case 
the  largest  error  is  generated  in  the  700-  to 
900-mb  region  and  the  error  reverses  sign  above 
the  100-mb  level. 


Subtracting  Eq.  (6)  from  (1),  after  using  (7)  to 
evaluate  the  integrals  and  rearranging  terms, 
yields  the  temperature  error 

* (At  - Ft  ♦ Cj  - Dt)/Eif  (9) 


?4o  260 

U r.*i ; i.atijpe  (®k) 


Figure  1.  Tmperatur e profile  used  for  oaloulationc . 
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Figure  4.  Retrieved  temperature  errors  for 
various  atmospheric  CO 2 deviation  profiles  of  the 
family  shotm  in  the  inset . a represents  the 
maximum  deviation. 


h +2  43 


cnroR(c:0 


TEMPERATURE  ERROR  (°K; 

Figure  6.  Retrieved  temperature  errors  for 
various  atmospheric  CC2  deviation  profiles  rf  the 
family  shown  in  the  inset,  n represents  the 
maximum  deviation. 

A.  DISCUSSION 

Fluctuations  in  the  CO^  mixing  ratio  below 
the  tropopause  have  the  greatest  effect  on  low 
altitude  temperature  errors  This  follows  from 
the  fact  that  most  of  the  spheric  CO2  is  be- 
low the  tropopause.  The  temperature  errors  at 
high  altitudes  caused  by  low  altitude  CO^  fluc- 
tuations (Figure  A)  are  the  result  of  the  smooth- 
ing technique  used  in  the  retrieval  program: 
each  of  the  scaled  constants  is  obtained  from  a 
weighting  of  all  of  the  radiometer  channel  resid- 
uals (Chahine,  1971).  The  net  result  is  that  the 
computed  temperature  of  any  altitude  Is  affected, 
to  varying  extents,  by  the  CO2  error  at  any  other 
altitude . 

As  an  indication  of  the  importance  of  de- 
generated errors,  temperature  errors  resulting 
from  the  use  of  averaged  0^  profiles  were  calcu- 
lated. Temperatures  were  obtained  from  the  re- 
trieval program  in  which  three  different  ozone 
transmit tances  were  used.  These  transmi ttances 
were  for  typical  equatorial,  polar,  and  averaged 
ozone  profiles.  The  differences  in  retrieved 
temperatures  using  the  average  profile  and  that 
using  the  equatorial  and  the  polar  profiles  pro- 
duced an  rms  temperature  error  of  0. 3°K  when  the 
entire  profile  was  averaged.  The  average  error 
below  500  mb  was  nominally  0.8°K,  with  a maximum 
error  of  about  1*K  near  ground  level.  As  can  be 
seen  from  Figure  3,  a 5%  CO2  variation  will  pro- 
duce a larger  temperature  error. 


• . * •m  en* urn  errors  for 

• » le  %;  it:  ~,n  profiles  of  the 


To  be  of  most  value,  statistical  expected 
temperature  errors  due  to  CO2  fluctuations  are 
needed.  This  can  not  be  done  at  present  because 
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3.  RETRIEVAL  METHOD 


u O' 

A 0*> 


' r,}2 
C02 


TCMPFRATURE  TRuOR  (*K  ) 

Figure  2.  Temperature  errors  determined  from 
Equation  (10)  for  various  uniform,  atmospheric 
C0£  concentrations.  CO o concentrations  are 
expressed  as  percent  of  320  ppm. 


These  results  can  he  easily  understood  if 
Eq.  (9)  Is  simplified.  For  small  a. , 


JPnP 


and 


AHnP,  * AlnPfl. 


(10). 


Fjcpandlng  B[T(Pfl)l  in  terms  of  B(T(P  )],  one  can 
obtain  a lowest  order  approximation  tor  the  tem- 
perature error: 


The  second  method  utilized  synthetic  data  in 
a temperature  retrieval  program  for  the  SIRS-B 
satellite  temperature  sounding  technique  (Duncan, 
1973).  The  transmi ttances  were  modified  for 
arbitrary  a^(P),  and  the  errors  between  the  re- 
trieved temperature  (for  the  modified  CO2  content) 
were  subtracted  from  the  standard.  The  results 
of  this  technique  are  shown  in  Figure  3 and  may 
be  compared  to  the  results  of  Eq . (10)  (Figure  2). 
Uniform  of  t5 % and  '10%  was  used  in  the  re- 
trieval program.  As  can  he  seen,  the  agreement 
with  the  closed  equation  solution  is  generally 
good . 
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TEMPERATURE  ERROR  (°K) 

Figure  3.  Temperature  error  versus  altitude  fer 
different  uniform  atmospheric  COo  concentrations 
determined  by  the  retrieval  program.  The  C0« 
concentrations  are  expressed  as  percent  of 
320  ppm. 
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To  a zero**1  order  approximation,  the  temper- 
ature error  is  simply  the  magnitude  of  the  local 
temperature  versus  pressure  slope  times  the 
pressure  difference  between  the  centroids  of  the 
standard  and  modified  weighting  functions.  The 
temperature  error  simply  results  from  the  fact 
that  the  sensed  temperature  is  for  a different 
pressure  than  anticipated.  In  altitude  regions 
with  large  magnitude  temperature  slopes,  one 
expects  large  errors;  where  small  slopes  exist, 
one  expects  proportionally  smaller  errors.  The 
error  changes  sign  with  the  change  of  slope. 


The  expected  temperature  errors  for  various 
CC>2  modification  profiles  are  shown  in  Figures  4 
through  6.  On  each  graph  the  CO2  variation  is 
shown.  The  parameter  u on  the  temperature  error 
curves  represents  the  maximum  CO,,  deviation  for 
the  particular  family  of  CO,,  deviation  profiles. 
For  example,  in  Figure  4,  a represents  the  CO2 
change  in  the  800-  to  1000-mb  region.  In  the 
range  from  R00  mb  to  100  mb,  n reduces  linearly 
from  its  value  at  800  mb  to  0 at  100  mb.  The  AT 
curves  are  for  the  particular  a chosen.  In  the 
other  figures,  a and  the  graphs  are  similarly 
re  lat  ecf. 


It  should  be  noted  that  Eq.  (10)  may  be 
utilized  to  examine  the  temperature  error  re- 
sulting from  any  absorber  perturbation.  One 
need  only  use  the  correct  combination  of  indi- 
vidual t ransmi t tances  to  correctly  specify  the 
total  transmittance. 
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the  amplitude  of  probable  CC^  variations  is  not 
known.  The  use  of  the  variations  in  monthly 
averages  or  the  variances  measured  at  ground 
level  seem  equally  invalid.  One  can,  however, 
compare  a 5%  C0o  fluctuation  temperature  error 
to  that  caused  6y  the  radiometer  random  error. 
The  reported  random  error  of  0.25  erg  sec”*  cm”^ 
sec”*  cm  (Wark  and  Hilleary,  1969),  which  when 
averaged  over  a given  locality  is  substantially 
less  (Fritz,  1970),  produces  direct  rms  tempera- 
ture errors  of  the  order  of  0.5°K.  It  is  clear 
from  the  figures  that  a 5%  CO2  uncertainty  has 
a larger  associated  error  in  the  lower  atmos- 
phere. It  might  be  noted  that  if  one  used  the 
lower  limit  approximation,  ij  = V1  + a^, 
the  basic  conclusions  are  unchanged. 
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Based  upon  these  results,  one  can  speculate 
that  a contributor  to  the  difficulties  in  equat- 
ing experimental  and  theoretical  15-pm  atmos- 
pheric transmit tances  can  be  a variation  in  the 
CO^  content.  This  is  not  to  say  that  aerosols 
or  computational  techniques  are  not  a problem, 
but  that  one  should  consider  the  possibility  of 
some  uncertainty  in  CC^.  One  can  equally  specu- 
late that  some  of  the  problems  in  temperature 
retrievals  rest  on  the  same  cornerstone. 
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where  Uj  is  the  mean  value  of  a(P)  over  the  layer. 

Sample  calculations  indicate  th.it  the  rela- 
tionship in  Eq . (A2)  is  a reasonable  approxima- 
tion for  the  15-um  transrai Ctances . Expressing 
the  modified  total  transmittance  In  terms  of  Eq. 
(A2) , 

- n At*.  ••  H At  : {All 


This  allows  the  modified  weighting  functions  to 
be  evaluated: 

(A4) 


31  i / 1 . -r  \ (Ho'  )3At 

’ (n  Vxr.  - — i- 

3lnP  P.p,  pklpj  pkipj  flUk  >mr  p-p 

J 


dtnAt. 

t*(P)U  (l«k)  ).  <A5) 


APPENDIX  2 


Homogeneous  Transmittance  Program 


78 


dimension  wi7),  r ( 3 2 s ) , g n u i 3 o o o » , s ( 3 o o o j , *LpHA(3noOj,  £ o p ( 3 o o o i 
DIMENSION  mol  (3000).  C A Y l ( 7 I , OpD ( 6000 ) , ENU ( * 000 ) , TRaNS(IOOO) 
DIMENSION  SUM  1 ( 7 ) . CS2 ( 7 ) 

DIMENSION  XNU I 325 ), WX ( | 00 » . ESUM(IOO) 

dimension  ww(inn) 

DIMENSION  QV ( 7 , 7 ) » TqV < 7 I , QVF  I 7 ) 

DOUBLE  PRECISION  CAYl.SUMl 

c this  program  generates  a transmittance  spfctrUm  with  Output  results 
C PRINTED  EVERY  deLV  WAVE  NUMBERS  BETWEEN  The  INITIAL  FHfQUENCY, 

c vi,  and  the  final  frequency  , V2.  calculations  are  perFormeo  for 
C A UNIFORM,  CONSTANT  PRESSURE.  CONSTANT  TEMPERATURE  PATh  CONTAINING 

c any  or  all  of  the  molecular  species  described  in  this  report 
c in  arbitrary  amounts,  molecular  abundances  must  be  specified 

C IN  THE  UNITS  ( M0LECULES/CM2 ) . MONOCHROMATIC  CALCULATIONS  ARE 

C MADE  At  FREQUENCY  I NTEr VALS , DV  , AND  A TRIANGULAR  slit  FUNCTION 

C OF  HALF-W I DTH , A , IS  CONVOLVED  WITH  THE  MONOCHROMATIC  RESULTS 

C 

c 

c 

DEPTH  = 0 • 00  1 
P I =3 . 1 4 1 59 
DO  8 6 *)  M = 1 , 7 

868  READ  < 5 , 862  ) ( Qy ( M , I ) , I ■ | , 7 ) 

862  FORMAT ( 7F 1 0. 5 ) 

TQV ( 11=175. 

DO  863  1=2,7 

863  T Q V ( I ) = T Q V ( 1-1  )+25. 

TQV ( 6 ) =296  . 

I R UN  = 0 

l 0 I S CONTINUE 
_ ' 1RUN=IRUN*1 

I EO  F = 0 
K E N D = 0 

. SUM-0.0 

I V=  I 

READ ( 5 , 305  ) I R .kfNDF 
DO  306  I = 1 , I R 

3 0 6 RE  AD ( 5 , 307 ) CODE  , XNU ( I ) , 9 ( I ) 

305  FORMAT (215)  ( 

307  FORMAT(F10.8,Fio.  1 ,F10.S) 

C 

C READ  INPUT  PARAMETERS  ( P = PRE5SURE  ) , ( T = T E M P E R A T U R E ) . 

C W ( l ) =H  2 0 , W(2)=c02,W(3)»03,W(8)=n2o,W(5)=c0,w(6)sCHH,W(71«O2. 

C VI  AND  V 2 ARE  FREQUENCY  LIMITS  FOR  WHICH  OUTPUT  RESULTS  ARE  REQUIRED. 

C DV  IS  MONOCHROMATIC  FREQUENCY  INCREMENT. 

c bound  is  the  frequency  prom  any  line  center  beyond  which  the  line 

c will  be  neglected. 

c A IS  The  HALF-WIDTH  OF  a TRIANGULAR  SLIT  FUNCTION, 

C DELV  is  FREQUENCY  INCREMENT  OF  CONVOLVED  OUTPUT  T R A N S M j T T A NC E 

C RESULTS. 

C 

READIS, asm, V2.0V, BOUND. A, DELV 

PRINT  87, VI  ,V2,oV. BOUND,  A. DELV 
V l 0-V I 


C IF( A*2/DV*1 .GT.3000)  CALCULATION  CannOT  Be  DONE 

C IF  THERE  ARE  HnflE  THAN  ->nOO  LIwEr  nnO  FBnH  TAPE  rN  A rREQUENCY  RANGE 


r 


i • 


• h 


/• 
)• 
?• 
3 • 
t • 
2. 
3 • 
4* 
5* 
S • 


1 1 


304 

700 


OF  2 ( A ♦ BOUND  1 CALCULATION  cannot  Be  done 
VBOT-V  | -A-BOUNO 
VT0P-V2* A^BOUNO 
1 = 1 
ILL=  1 

REA0(2,9S)GNU(l)  ,S<I),ALPHA(I)»EDp{I>iIDaT,IS0T.M0L<1> 
I F ( GNU ( II  .L T . VfiOT > GO  To  t 
IF (GNU< I ) ,GT. VTOP >G0  To  11 

I F ( MOL  C n . NF  , 3 ) GO  TO  J 
I F ( 1 . GT . 2999  ) Go  TO  1 1 
1*1*1 

GO  TO  1 

II  = 1 

READ(5,77)P,T 
PRINT  79,  P,T 
REA0(S,8! I <W(Mj ,M=I ,7) 

PRINT  83 

PRINT  0 I , C VK  ( M ) ,M=  I , 7 ) 

DO  304  M»1 ,7 

V*  < M 1 = N(  M 1 »2 , 69e*  1 9 

PRINT  700  , (VMM)  ,M=1  ,/) 


F0RMAT(|4H  W IN  M0L/CM2.7(2X,Et2.6)  > 

SET  NX  VALUES  (MULTIPLES  OF  THE  N<I>  VALUES) 
NX ( I ) = 1 .0 
DO  1000  K = 2 , I 00 

. NX(K)=NX(K-I  )*S.0-  - 

1000  CONTINUE 

min*o 

DO  870  M = 2 , 7 
M l N = M I N* 1 


7*  . 
8* 

9* 

0. 

1 • 

2* 

3 • 

4*  _ 

5. 

6 • 

7* 

8* 


IF(T.LE.TQV(M) )gO  To  071 

870  CONTINUE  - 

871  DQV=(T-TQV(MINn/(TQV(MINtl)-TQV(MiN)) 

00  869  M = 1 , 7 

869  QVF(M)-(<}V(M,Min)*(dQV*(qV(M,MIN*1)-QV(M,MJN)>))/QV(Mi6) 
KEND=KEND+1 
DO  1003  K=  1 , 1 Oo 
1003  FSUM(K)=0.0 
FSUM1=0.0 

PRINT  97,  VBOT  , vTOP ,GNu ( I 1 1 , I 1 
15=1 

V2P  = GNU  (111  -BOijnO-A 


9* 

0. 

; 1 • 
1 2 • 
')3« 


1- 

ji: 

■ f)7» 
3 8 • 
39. 
. 10. 
1 I* 

f12- 

r 


■ 


|C  ' 

C TAPE  has  BEEN  read  for  all  NECESSARY  lines  or  for  THE  MAXIMUM  NO. 

C OF  LINES  POSSIBLE  subject  TO  RECYCLING* 

C HALFNIDTHS  NILl  BE  SUPPLIED  BELON  when  they  DO  NOt  APPEAR 

C ON  TAPE. 

c 

00  15  I = I LL  , I 1 
M-MOL ( I ) 

IF  ( M . EQ • 1 1 GO  TO  IS 

IF  ( ALPHA ( I I .Gt.O.O)  Go  TO  13'  - 

IF  (M.EQ.2)  ALPHA ( I ) «0.o7 
IF  1M.E0.31  ALPHA!  I 1=0.  | 1 
IF  (M.E0.41  ALPHA ( I 1 =0 , 08 
IF  (M.E0.5)  ALphA ( I 1=0, q6 
IF  (M.EQ.6)  ALPHA ( I 1 =0 .055 


, ■ , 


80 


25. 

DO 

2 1 M=  1 ,7 

2 6. 

IF 

( M . EQ • 1 ) 

GO 

TO 

1 7 

27. 

I F 

I M . EQ . 2 ) 

GO 

TO 

1 9 

28. 

IF 

IM.EQ.3I 

GO 

TO 

1 7 

29. 

IF 

( M . EQ . 9 ) 

GO 

TO 

1 9 

30. 

I F 

(M.EQ.5) 

GO 

TO 

1 9 

3 1*  - 

IF 

( M.EQ. 6 ) 

GO 

TO 

1 7 

32. 

IF 

( M. EQ . 7 ) 

GO 

TO 

1 9 

33.  1 7 

CS2 ( M ) = ( t TO/T  ) 

• • 1 • 

5) 

3 R • 

GO 

TO  2 1 

IF  (M.E0.7)  ALPHA  (I  ) "0.  nM8 

IF  ( ALPHA!  I ) .LT.0.01  .OR.ALPHAI  I ) .Gr.  I .0)  ALPHA ( I > «0 . 06 

CONTINUE 
I S°  1 

P0=  l 013.00 
T0-296.00 

CSI-(T0-T)/(T0*t*0.69H6) 

ROTATIONAL  PARTITION  FUNCTION  IS  DEFINED  BELO* 


CS2(M)»T0/T 
CONTI NUE 

C A = < (T0/T)**0.5)*tP/P0) 

TEMPERATURE  DEPENDENCE  OF  all  LINE  INTENSITIES  COMPUTED  HERE. 

V=V 1 -a 
DO  27  M=  1 , 7 
CAY  I (M  ) =0.0 
SUM  1 <M)=0.0 

DETERMINE  INDICES  < 1 5 AND  16)  INDICATING  WHICH  SPECTRA^  LINES 
ARE  TO  BE  USED  jN  THE  CALCULATION  AT  FREQUENCY  V. 

DO  33  1 = 15, I 1 

,IF  ( V-BOUND-GNui  I ) ) 29,29,33 

I So  I . . .... 

GO  TO  3 5 

CONTINUE 

I 5 ■ l l 

GO  TO  H 9 

DO  39  Jo  I S , I 1 

IF  ( V + BOUND-GNu ( J ) I 37,37.39 
I 6 = J- 1 
GO  TO  M3 

CONTINUE  < ' i 

16=11 

COMPUTE  THE  optical  DEPTH  AND  TRANSMITTANCE  at  FREQUENCY  v. 

DO  M 5 1 = 15,16 

M=MOL ( I ) 

Z=ABS I V-GNU ( I ) ) 

X1«S(  1 )*CS2(M)*EXP(_E0P(  I ) • C S I ) • Q V F ( M ) 

X2-ALPHA I I ) *C A 

SUM l CM)»Xl*X2/(Z**2*X2*«2> 


,i . ■ - 


CAYl  (M)  * C A Y I ( M ) ♦SUM  1 <M) 

4 5 CONTINUE 
C A Y * 0 • 0 
DO  4 7 M=  1 ,7 

4 7 CAY  = CAY*CAYI  (M) ,W(M) 

OPD ( I V ) *CAY*0. 3 1 83 
GO  TO  51 

49  OPD I I V ) *0.0 
5 1 OPD ( I V ) =EXP ( -Opo ( IV)) 

IF ( V.GT.XNUI  I R ) )G0  TO  330 

335  IFIV.GE.XNUI  IF)  . A ND . V . L E . * NU < I F M 1 1G0  TO  337 
IF  (V.GT.XNUI  I F ♦ j ) )G0  To  336 

IFIV.LT.XNUI  IF)  ) 1F-IF-1 
GO  TO  335 

336  I F * I F * 1 
GO  TO  335 

337  DRF»R(IF)*<V-XNU<m>*(R(IF*l)-R(lFJ)/tXNU<IFM>-XNUllF>) 

DO  1001  K'l,10o 

1 0 o I FSUM ( K ) -FSUM ( K ) + D R F • I 0 P D ( I V ) • • W X < K ) ) 

FSUM I«FSUH  l ♦DRF 
GO  TO  339 

338  WR I TE < 6 , 340  ) V 

340  FORMAT ( 5X ,9HN0  A D D , V = i E 1 4 • 6 ) 

339  CONTINUE 

IF  | ( V + DV)  .GT.V2P)  GO  TO  53 

IF  ( V . GE  . V 2 + A ) GO  TO  53  — 

IF  ( I V.GE.6000)  GO  TO  53 
I V * I V ♦ I 
V * V ♦ D V 
GO  TO  25 


AT  THIS  POINT,  CYCLE  BACK  TO  STATEMENT  25  AND  COMpUTE  THE 
MONOCHROMATIC  T R ANS M I T T A NC E . A T V*Dv,ETC. 

IF  STATEMENT  53  IS  REACHED  ,ALL  POSSIBLE  MONOCHROMATIC  TRANSMITTANCE 
VALUES  HAVE  BEEN  COMPUTED.  AND  ThE  slIT  FUNCTION  CONVOLUTION  WILL 
NOW  BE  PERFORMED  IN  LOQo  57 


F R E Q = V 1 

PRINT  101,  I V , v , V 2P 
WRITE  OUT  TRANSMISSIONS 
WRITE(6,201  ) ( Opo ( L 2 ) ,L2» 
FORMAT!  1 OF  1 2 . 8 ) 
FINAL=V1*6000.*oV-A-DElV 


MONOCHROMAT I C 
I V ) 


JFNU=  l 

L=0ELV/DV>0.01 
I A = 1 

SUM»0 . 0 

SUM2=0.0 

DO  57  I = I A , I V 

SUM  = SUM.*<A-ABS(v-FRpQ)  ).0PD!  I ) 
SUM2»SUM2* ( A-Aps ( V-FREQ  ) ) 

V« V*0 V 

IF  (V-(FREQ*A>)  57,59,59 
CONT I NUE 

59  TRANS C JFNU ) *SUm/SUM2 
FNU ( JFNU 1 ■FREO 


— 

I~—  

c 

_ c 

c 

1 

c 

a 

1 - 

c 

C , 

\— 

5 3-— 

c 

I ...... 

L 

2 0 1 

V 

- 


ZB* 

IF  (FREQ.GT.V2)  GO  tO  61 

29. 

IF  (FRE0.GT.V2pj  GO  TO  61 

30. 

IF  (FREQ. GE. FINAL)  GO  TO  <»  1 

3 l • 

freq=freo*delv 

32. 

IF  (JFNU.GE.10oo)  Go  To  6l 

33. 

JFNU" JFNU* 1 

39. 

I A a I A * L 

3 S • 

V"FREQ-A 

3 A. 

SUM  = 0 .0 

37. 

GO  TO  55 

38* 

C 

39. 

c 

CONVOLVED  TRANsmITTaNCE  RESULTS  ARE  NOW  PRINTED  OuT. 

HO. 

c 

M 1 . 

61 

PRINT  103,  JFNU 

H 2 . 

PRINT  109 

H 3 • 

PRINT  105,  (FNU( J)  ,TRANS( J)  , J-l  .JFNU) 

HH. 

IF  (FREQ.GE.V2)  GO  TO  75 

H 5 . 

IF  (FREQ.GT.V2p)  GO  TO  67 

H 6 • 

IF  ( JFNU.GE. 1000)  Go  To  65 

H7. 

IF  (FREQ. GE. FINAL)  GO  TO  *3 

H 8 • 

GO  TO  75 

H 9 . 

63 

V 1 =F INAL.DELV 

50. 

15=1 

3 1 • 

I V * 1 

32. 

JFNU" 1 

53. 

V " V 1 - A 

39. 

GO  TO  25 

35  . 

65 

I A " I A * L 

;6. 

C 

57. 

C 

IF  STATEMENT  6s  IS  REACHED,  ADDITIONAL  MONOCHROMATIC  CALCULATIONS 

58. 

c 

are  required  TO  SATISFY  the  TOTAL  FREQUENCY  RANGE  OVER  which 

39. 

c 

CONVOLVED  RESULTS  ARE  REQUIRED. 

so. 

c 

3 1 • 

JF  NU= 1 

S 2 . * 

v=frfq-a 

S3 

s9 

55 

3 t> 

S7 

S« 

,9 

I '0 
r I 
i '2 

j 7 3 


GO  TO  55 
67  IV=1 
C 

C IF  STATEMENT  67  IS  REACHED.  THE  0*TA  FROM  THE  DATa  TAPe  WILL  Be 

c reorganized  and  the  tape  will  be  Read  again. 
c 

J F N U = 1 
V 1 a F R E 0 

VBOT=V 1 -A-BOUNd 
DO  69  IN  " I , I I 

IF  (GNU(  IN)  .GT.vBOT)  Go  TO  71 
69  CONTINUE 


I N = I 1 
7 1 I J 3 I N 
L - 1 

00  -7  3 I = l J , I 1 
GNU ( L ) ■ GNU ( I ) 

S (L ) «S ( I ) 

ALPHA ( L ) "ALPHA ( i ) 
EDP ( L ) "EDP ( 1 ) 
MOLfL)"MOLl I ) 

L"L*  1 


— — 


I 


I 


» 

i 

) 

1 

7 

) 

3 

3 


1»L 
ILL-L 
GO  TO  1 
75  CONTINUE 

DO  1004  K= I , 1 Oo 
F SUM ( K ) = FSUM ( K ) /FSUm  1 
WU  = W ( 3 ) *WX  < K > /2 . 69E  1 9 

WW  ( K ) =WU  • 

WRITE<6,341)C0D£»P.T.WU,FSUM(K> 

1 Oo  4 CONTINUE 

»*R  I TE  < 7 , 2002  ) ( WW(K  ) ,K=  1 . J 00  ) 

2002  FORMAT ( 8F I 0 .5  ) 

WRITE(7,2003)(FSUM(K),K=1.100) 

2003  FORMAT ( 8F  10, 7 1 

341  FORMAT { 30H  03  UNIFORM  TRANS.  FOR  I nST. ,F I 0 . R , I 3h  AT  PRESSURE  ,F 8 

1 • 3 , 1 0 H AND  TEMP,  ,F8,3, 1 ?H  03  A TmCm . • 1 P?  . 5 , l 1 H TR A N S M I T . • . E 1 5 . 8 » 

IF (KEND.GE.KENDF)GO  TO  202 
1 « 1 1 
V 1 = V I 0 

V = v 10  - - 

I V = 1 
I L L * 1 
I E OF  = 0 
SUMaO.O 

VBOT=V 1 -A-BOUNq 

VT0P'=V2  + A + B0UNd  — - - - 

GO  TO  11 
202  CONTINUE 
REWIND  2 

I F ( IRUN.LT. 6)Go  TO  1 0 1 S 
CALL  EXIT 


s. c 


;• 

77 

format 

(E12.5.F7.2) 

3 • 

8 1 

format 

I 7 E 1 0 . 3 ) 

?• 

85 

FORMAT 

( 6F 1 0 . 3 ) 

j* 

79 

format 

('  PRESSURE  =*  ,e 1 ?•&, * TEMPERATURE  ■ • , F 7 , 2 ) 

1 • 

83 

format 

(3X,'WATER,.6X,*c02»t6X.,OZ0ME,»7X,'N2O«,7x,'CO.,8X,fCH4',7 

2* 

1 X , * 0 2 ' , 

4 X 1 

3* 

87 

FORMAT 

(•  vi  =»,fio,3.*v2  = • , f 1 0 , 3 , • ov- * , f i 0. 3 , • Bound  »*,fio.3,*a 

4 • 

l3’  ,F10. 

3,’DELV  =•  1 F 1 0 • 3 ) 

5* 

89 

format 

<•  PARITY  ERROR  ENCOUNTERED  AT**FI2.3) 

6 • 

1 9 1 

format 

(*  END  OF  FILE  ENCOUNTERED »,  I 5 > 

7. 

97 

FORMAT 

(*  VBOT  = * ,F l 2.3  , * VTOP  ■* ,Fi2.3,»GNU  »,,Fl2.3.'ll  • • i S > 

3 • 

103 

FORMAT 

( • JFNU  c • , I 5 ) - 

9* 

1 04 

FORMAT 

I 5 ( • FREQUENCY  TRANS.  •)) 

( 

0 • 

93 

FORMAT 

(110) 

1 • 

95 

FORMAT 

(F10.3,E10.3,F5.3|F10.3,3SX,I3,1!1,I3) 

2* 

l 0 1 

FORMAT 

( I 5 , 2 F 1 0 , R 1 

'3* 

1 05 

FORMAT 

(5IF10.3.E12.5)  ) 

4* 

STOP 

5* 

END 

1 1 

1 

lor 

compilation: 

NO  DIAGNOSTICS* 

i woo;  *c H r. c k • n o n i r . w , » 

W A S G i F T • • T , I Ml  4 . H ;;  F L 3 1 7 1 7 

l»'  A S H , T 2 • , F ///  SOU 
fcCOP  Y , m T . , 2 
iS  F K L L T 

13  MS(.,  .V  COlUlU  f-  O 0,'UCt  PAN  8 2U0 
® M 5 Vj  , »J  P L S E SeNi)  C 0 N r R f A p E T 0 BLDG  i A ? 3 
'»  a s g , r q , t , C o 1 o i o 

® F 0 R , I S • U N I F 

0 | MENS  I ON  «.•  t 7 ) , ;,NU  I 3000  ) , S ( 3000  ) , ALPH  { 3 iQ'>  ) , F.  0 P ( 3 0 0 , . ) 

0 ! M E N S 1 0 N U P 0 t 3 OH  111  , F N 1 1 ( 1 II 1 1 0 ) , 1 K A N S ( 1 n U 0 ) , X L ( 3 0 U ij  ) , C 5 2 ( 7 ) 

• > I M t N s l 0 N w V ( 7 , 7 ) , r 10  V < 7 I , u V K ( 7 ) 

U | Mt.  N , ( ON  SYIIH  ( ,,  ) 

DOUBLE  PRECISION  CAY,  SUM  l,  V 10,  V20 

DOUBLE  P K fc  C I S 1 0 M i<  V X , v X , I-  R E 0 X , 0 E L V X , A , , X r>  , V 1 X , F 1 U A L X 

c 

c 

C This  PROGRAM  cC.'Ei  at.  s A TPANSO  I TTANCr  SprCTRUM  „ 1 T H nUTPUT  «r~,'JLTS 

C *’  F;  INTEL)  L V t !<  Y Of.U/  -/A  V r NUMBERS  0 F T A t N THE  INITIAL  FrEQUENdY, 

C V|',  AND  The  F(NaE  FK  E )UEnC  Y , V 2 . CALC;  LA|  IONS  A R f.  PERFORMED  FOP 

C A UN  I F OKM  ,, CONSTANT  PRESSURE,  CONSTANT  TEMPERATURE  path  CONTAINING 

C Amy  OK  ALL  OF  The  MOLECULAR  species  0 1 S C r I B E 0 IN  THIS  REPORT 

C In  ARBITRARY  AMOUNTS.  MOLECULAR  ABUNDANCES  MUST  BE  SpECiFIeO 

C In  THE  UNITS  (MoLECUl.ES/C  M2)  • MONOCHROMATIC  CALCULATIONS  ArE 

c made  at  frequency  intervals, ov,  and  averaged  over  a frequency 

C I nTER  VAL  = 2A  .LOReNTZ  LInL  SHAPE  IS  USEp. 

c 

C This  PROGRAM  USES  Tape  This  ( 3 1 7 I 7 ) Oe  FasTRAN  file  «200*HCOti, 

c For  input  data* 

c 

c The  output  is  to  tape,  specifies  pressure,  temperature,  gas 
c anq  averaged  tRansmittances,  starting  wavenumber  and  interual.dElv 
c and  number  of  data  points. 

c 
c 

c S 0 = 1 . c 

C A 0 = I .0 

1 E 0 F = C 
OePTH  = C].001 
p I = 3 . 1 4 I S9 
S u M = 0 . 0 

I V = 1 

Rewind  a 

c Read  number  of  layers 

KEA0(B,8  7 2)  ILAY 
8 7 2 F 0 R M A T ( l 5 ) 

I L A = I 

* R I r E ( 0 , 1 0 7 ) 

Wtfl  T E ( B , 0 7 3 ) 

873  Format {♦NUMBER  of  LAYERS  for  calculation  FOLLOWS*) 
r I t E ( H , 8 7 2 ) l L A Y 

c 

c Read  partition  fn.oata 

c 

DO  864  Me  | , 7 

864  Read  <S,&62><QV(M,I),  1 = 1,7) 


062  r OKU  A T ( 7 M U . 5 ) 

T R V ( 1 ) = 1 7 b • 
l>0  H6J  1=2,7 
0 6 3 T i)  v ( I ) = T 0 V ( I - 1 ) * 2 b • 

T u V ( <>  ) - 2 V 6 . 

C 

C K<  AO  I Ni’llT  PARAMETERS  ( I’ - ),  ( r ~ T | M f ’ L i(  A T U K t I i 

C «'•  ( i 1 - H 2 0 , .v  ( 2 ) = c <j  2 , VV  ( j I = o 3 , W < H 1 - N 2 0 , IV  ( b ) = C 0 » 77  ( A ) = C H 9 , rt  ( 7 ) = 0 ? . 

C V|  AN(>  V 2 A PL  F,;t.wULMCY  LIMITS  f HR  vViirCH  OUTPUT  RESULTS  ARE  REQUIRE,)* 

C DV  IS  MONuCmRoMaT  l C FR’F.QULNCY  INCREMENT  » 

C I'OUMO  IS  THL  FREQUENCY  I ROM  AMY  LIME  (ENTER  REYOilD  <<  H I C H THE  LlNf. 

c will  he  nfglfc Tpo . 

C a IS  The  HALF- WIDTH  OF  A TRIANGULAR  S i IT  FUNCTION. 

C o I L V IS  FREQUENCY  INCREMENT  0 F CONVOLvLO  OUTPUT  TRANSMITTANCE 

c Results*  ' 

c 

C For  C 0 2 FILTERS  NCAA  OR  PAP,  Vl=605.0,VZ=767.0,DV=.00l,DELV*.Ul 

C buUND  = 1 5 . (J  , A = . OpS 

c 

Read ( s ,096 ) mw  ,sn  in 
89  6 format  ( i s , l i £'. 3 » 

REA0«S,e5)Vl,v2,DV,B0UND,A,0ELV 
0 E L V X = 1 .00-2 
OeL  V = OELVX 
A x = b . OD- 3 
A = A X 

0 v X = l .00-3 

0 v - 0 V X 
v 1 X = V 1 

v 1 0 = V 1 X 

V 2 0 = V 2 

Print  07,vi ,V2iOV,bound,a,delv 

C l FI A*2/0V+ 1 . GT • 3000  I CALCULATION  C A N N n T BE  DONE 

C IF  THERE  ARE  MORE  THAN  3000  LINE'S  RE  A n FRQM  TAPE  IN  A FREQUENCY  RANgE 

C . Of  2 < A + 8 0 U N D ) CALCULATION  CANNOT  BE  DpNE 

V B 0 T = V 1 -A-bOUljD 

V T 0 P = V 2 * A * b 0 U N D 

1 = 1 

I LL  = 1 
R E V.  1 U u 2 
I RE A0  = 0 

1 READ12,9S)GNU|I),SII),ALPHA(I),E0PII).IDAT.IS0T,M0 
I F (GNU ( I ) ,LT . VBOT ) GO  TO  1 
lF(GNU(!).GT.VT0P)G0  TO  II 
1 F ( M 0 . N E • N V.  > G 0 TO  1 
I F ( S C I I .LT • SM I N ) GO  To  1 

lF(Et»H(I),LT.&Oo.)GO  TO  6 50  1 

I F < S < I ) • GT • 1 • OE - 2 3 > GO  TO  650 

IF(EOP(I).LT.IOOO.O.ANO.S(I).GT.I.OE-?H)GO  TO  650 
lF(GMU(I).GT.72o..ANO.S(I).GT.I.OE-25)GO  TO  850 
Go  TO  1 
850  Continue 

1 p ( I . G T • 2 9 9 9 ) (,  0 TO  11 
I = I ♦ 1 
<’0  TO  1 
II  I l ■ I 


87 


— 


4 


211 ‘I 


c 

c 

c 

c 

c 

c 


1 3 
1 S 
C 
C 

876 

780 

777 


789 


700 

C 

C 


8 70 

07  i 

869 

8 7 8 


N 0 • N ft 
l nil  1 I Mill 

I’kINI  77,  V It  0 T t V Till’  , f,fllJ  ( I 1 ) , I l 
10=1 

V ?(•  = (, Nl)  ( 1 I ) - HoUnD-  A 
If  ( V 2 P . L T , V 2 ) (,  0 TO  8 7 0 

T A P L HAS  BEEN  K f.  A 0 FOR  « L L NECESSARY  | 1 N y 5 OR  FOR  THE  M A x I M u M MO. 

Of  lines  possible  subject  to  recycling* 

half  widths  /.ill  he  supplied  below  when  They  do  not  appear 
On  tape* 

Oo  is  I - I L L , I 1 
N - N tl 

If  ( 1 1 . F Ci  . 1 ) fi  0 TO  IS 
If  ( alpha  ( i i . CiT  . u.n  > c;o  to  i 3 
If  ( M.  Ed  . 2 ) ALPHA  ( I ) =0  . l>7 
If  ( M . E G . 3 1 A L P H A l I ) = 0 . 1 | 

If  ( N . F 0 . H ) A L P H A ( I ) = 0 . 0 8 
If  ( M . E 0 • S ) A L P H A ( I ) = 0 . U 6 
If  < M . E y • 6 > A L p H A I I ) = 0 . 0 s S 
If  ( H . L 0 • 7 1 A L P H A ( I 1=0.088 

If  ( alpha ( 1 ) .lT.0.01  .or. alpha ( I ) ,GT.  1 .0)  ALPHA < I ) «=0  • 06 
Continue 

Read  layer  data 
Conti nul 

Do  7 80  11=1,1000 

TrANS(M)=O.U 

R E A 0 ( S , 7 7 7 ) P I , T 1 ,P2  , TP. 

F0RHAT(2<e12.5,f7.2)  ) 

P= (Pi + p 2 ) / 2 . 

T=(Tl*T2>/2. 

print  7E9,pi,p2,p,ti,t2,t 

F0RMAT(2A,'Pl=»,F7.2,'P2=«,F7.2,'PAV=ttF7.2,'Tl=',F6.2,*T2»i,F6»2, 

* TAV • , F 6 . 2 I 
Op  = P 2 - P 1 

R E A 0 1 S , 8 1 ) { ft  ( M > , M = 1 , 7 ) 
print  a 3 

Print  a i , ( ( t t ) , m = i , 7 ) 
print  7 0 0 , < a < m ) , m = 1 , 7 ) 

Format cirh  w In  mol/cm2 , 7 < 2x  ,e  1 2. 6 ) ) 
modify  bound  for  pressure • ll . 3 . s *p 

B0UND=,01m*P 

If (BOUND*  LT.O.QOSIBOUNDsO.OOS 
M I N = U 

Do  8 7 0 M = 2 , 7 
M I N = M l N + 1 

lF(T.LE*TQV(M)UO  TO  871 
C ON  T i hUt 

D 'j  V = ( T - r :)  V ( M I N ) ) / ( T J V ( M I N ♦ 1 ) - T 0 V ( M I N ) , 

0 0 8 6 9 M = | , / 

0 V F t M ) a ( 0 V ( 1 1 , M I N ) + ( 0 W V • ( G V ( N , M I N ♦ 1 ) - Q v ( M . M I N ) ) ) ) / Q V ( M , 6 ) 

*<  W I T C ( 8 , 8 7 '*  ) 

F 0 K M A T ( • F 0 R L A Y E R NUMBER*) 


I 


I 


88 


fifcru 


AD-A069  041 


UNCLASSIFIED 


TEXAS  UNI V AT  EL  PASO  DEPT  OF  ELECTRICAL  ENGINEERING  F/G  4/1 

ATMOSPHERIC  TRANSMITTANCE  STUOY  WITH  THE  METEOROLOGICAL  SATELLI— ETC<U) 
OCT  76  R E BRUCE*  J H PIERLUISSI*  S K WEAVER  DAEA18-76-C-0019 


PR4-76-DC-30-PT-2 


NL 


END 


7 -79 


>v  K I T K ( H , *>  7 2 I I l.  A 
w K I T I.  ( M , I U / I 

107  J-'OKHA  r ( • MdHOfj  AH(  OU'J  PATH  OUTPUT  FOP  G.lSES  AS  FOLlOaS*  ) 
w win  ( n , *1  1 ) I U I U > , M = | , / ) 

Ar  I TL  Ml  ,»»!.»  V I , v,>  , UV  , HOUND  , A , OCL  v 

ftw  i ni  n . / 7 v ) p i ,p2  ,p  , 1 1 , 1 2 , r 

7 7 9 f o R M A T ( 6 ( 2 A , F 7 . 2 ) I 
X = ( V «>  V l ) / OIL  V + U • 1 
N X = ( A * I .0  ) 

R I If.  I R i 1 0 8 ) N X 

108  Format  < • number  of  data  points  =',I10) 

I 5 = l 

H 0 = l u 1 3 • n 0 
T o = 2 9 6 . 0 U 

Csl=(T0-T)/(T0*T*0.&9H6» 

c 

c Rotational  partition  function  is  defined  s r l o a 


Do 

2 1 M= 1 , 7 

IF 

( M . EQ . 1 ) 

GO 

TO 

1 7 

if 

( M . EQ • 2 ) 

GO 

TO 

1 9 

if 

( M . EQ • 3 ) 

GO 

TO 

1 7 

if 

( M . e y • M ) 

GO 

TO 

1 9 

if 

( M . E Q • S ) 

GO 

TO 

1 9 

if 

( M . E Q • b ) 

GO 

TO 

1 7 

I F 

( M . E Q . 7 1 

GO 

TO 

1 9 

Cs2(M)=((T0/T) 

* • I . 

5 ) 

Go 

TO  21 

19  Cs2(M)=T0/T 

2 1 CoNTl  / I U E 

C A - ( t T 0 / T ) * * 0 . 5 ) * { P / P 0 > 

CAX=CA/CAO 

Do  882  1=1,11 

082  ALPHA(I)=ALPHA(I>»CAX 
C A 0 = C A 

' SETUP  XL ( 1 ) 

0 0 ee3  i = i,ii 

883  XL(  1 )=$<  I )*CS2(mWD»QVF(N,v)*£XP(-E0P(  I i*Csi  I 

temperature  dependence  of  all  line  intensities  computed  here. 

vx=vix-ax 

v = vx 

2S  CaY=0.0 

determine  indices  (is  and  i&>  indicattng  ivh i cm  spectral  lines 
Are  TO  bE  USLO  in  CALCULATIONS  AT  FREQUENCY  V. 

DO  3 3 I = I S , I 1 

If  (V-F.OUND-GnIMI)J  29,29,33 
! 9 I S = I 

Of  TO  35 

13  Continue 

1 s = i i 

Oo  TO  99 
Do  3 9 j = l L , I I 


M M C M ,4 


2 9 C H A r 


S 


Ik  ( V ♦ t'OUNU-fiMU  I J J ) J 7, 3 /,3  V 
37  i 6 = J-  l 

t’O  r U M3 
3 V Conn  NOE 

I A - I I 

c 

C COMPUTt  THt  OPTICAL  DEPTH  AND  T R A N S M I t T A M C E AT  FREQUENCY  V. 

C t USES  L 0 K L N T Z LINE  ONLY 

C 

M3  0 0 Mb  I = I b i I 6 

i - A ‘J  S ( V - G N U ( 1 ) ) 

I F ( Z . G I . (3 C U N D ) 6 o TO  a 5 1 

St,iN|=Xl(I)*ALPHA(J)/(Z»Z*ALPHA(I)»ALPnA«J)J 
Go  T 0 l,  b 2 
d 5 i S U M 1 = 0 . 

65 2 CaY=CAY*SUM| 

ms  continue 

0 P D ( I V ) = C A Y • 0 . 3 1 a 3 • *.  ( N .v  ) • D P 
GO  TO  51 

M 9 OpQ(  I V ) = 0 , 0 
5 1 OpD ( 1 V ) = E X P I -OPD  I I V ) ) 

339  Continue 

lp  ((V+DV1.GT.V2P)  GO  TO  53 

1 F ( V . E . V 2 + A ) GO  TO  53 
If  (IV. GE. 3000)  GO  TO  53 

I V=  I V*  ! 

Vx=VX+DVX 

V = VX 

GO  TO  25 
C 

C At  THIS  POINT,  (.YCLE  BACK  TO  STATEMENT  2$  AND  COMPUTE  THE 

c monochromatic  transmittance  at  v*dv,Etc* 

c ip  statement  s3  is  Reached  ,all  possible  monochromatic  transmittance 
c values  have  been  computed,  and  the  slit  function  convolution  ,vill 

C . Nort  t)L  PERFORMED  IN  loop  57 

c 

5 3 F X E a = V 1 

F R E w x = V I X 

f H E U = F R E 0 X 

print  i o i , 1 v , v , v 2 p 

C ORITE  oor  TRANSMISSIONS,  MONOCHROMATIC 

20 l E ORMA  T ( l OF  I 2 . a 1 

Xn=3.UD*3«OVX-Ax-OELVX 
f I M A L X = V 1 X ♦ X 0 

final  =finalx 

Vx  = V 1 X-AX 

V r.  v X 

'•  j F N 0 = 1 

L = 0ELV/DV-m}.01 
I A = 1 

5 5 s U M = 0 . (j 

S 0 M 1 = 0 . 0 
Do  '3  7 I = I A , I V 
S 0 M = S U ! 1 * 0 P D ( I ) 

5 0 M l = 5 U M I ♦ 1 .0 
Vx»VX*DVX 


c 

c 

c 


c 

c 

c 

c 


5 7,5  7 


V = V X 

( V - ( F K I-  lj  4 A ) I 

C ON  r INul 

I K A II 5 t J F mj  ) = S u M / S u M I 
t NUIJCNUIstHij, 

If  ( F H L Q . fa  I . V 2 I vi  0 TO 
I r (M<tw.  fat. FINAL)  fa  0 
F R E >3  X = F R E Q X M)  E L V X 
T RCn=FkEUX 
If  ( J F M U . 6 E 
Jfno=jfnu+i 

1 A = I A +L 
Vx=FREr„X-AX 

v = vx 


6 l 
TO 


6 I 


lOOO)  GO  TO  61 


SijM=0 
GO  TO 


0 

5 5 


convolved  transmittance  results  are  Nnw  printed  out* 

6i  print  103,  jfnu 
print  i or 

print  ios,  ( fnu i j ), trams < j ), j= l , jfnu ) 

w R I T E ( a , 1 U 9 ) J F N u 
10V  FORMAT! IS) 

write (hi  trams 

106  F 0 H IT  A T ( 6 E 1 2 . 6 ) 

IF  (FREU.GE.V2)  go  to  7 s 
If  (JFnU.GE.IOOo)  go  to  65 
If  (FKEU.GE. FINAL)  Go  TO  63 
Go  TO  7 S 

63  Vix=FREUX*UELVX 
1 V 1 = V 1 X 
I 5=  1 
I V = 1 

Jfnu= i 
Vx  = v 1 X-AX 
v = vx 

Go  TO  25 
65  I A=  I A + l 

If  STATEMENT  65  IS  REACHED,  ADDITIONAL  MONOCHROMATIC  CALCULATIONS 
Are  R E (3  u I r e d to  SATISFY  thf.  total  FREQUENCY  RANGE  OVER  Which 
convolved  results  are  required. 


1 ■*, 


JF  NU=  1 

vx=freqx-ax 

V = VX 

go  to  ss 
75  continue 
jfnu  = u 

"H  * TL  I <1  , I ov  ) jFnU 
I L A = 1 L A ♦ 1 
J F N U = 1 
I 5=  I 


■MiHmps*, 


R t A 0 ( H , I 1 Cl  ) < 5 Y I I LM  I ) , I = I , 3 ) 

K E A J « M , » 7 2 ) I L A 
«E  AiH  b , 1 10  ) ( SYMh  ( I ) , I = 1 , H 1 
110  F o K M A T ( 7 A fr , A 2 ) 

WR l TE (6  , I 101  < 5yMH  (11.1=1.3) 

>V  R I T L‘  < 6 . b J ) 

R £ A 0 ( 8 , B 1 ) l .V  ( M 1 , M = 1 , 7 1 
WRITEl6.bl)(V\(M),M=l,7) 

REAOtb  i 6b  ) V 1 ,V2*0V,liOUN0,A,0ELV 
t>H  I TE  < A • «7  1 V 1 , V2 ,0V  , oOiJNO  , A , DEL  V 
l<  E A 0 < 8 , 7 7 V IP  t , P / , P , T l , T 2 , T 
V»  R 1 T E ( 6 , / B V 1 0 l , P 2 , P , T l , T 2 , T 
READ  ( 8 , 1 Ci  « 1 N x 
A R I T E ( A , 1 0 8 > N X 

500  Read  « n , i o v » j f n u 
AR 1 T E ( 6 » 1 0 9 ) J F N U 
lF(JFNU»Efii.01G0  TO  50? 
lF(JFNU.LE»NX)Go  TO  501 
WRITE(6,112) 

503  continue 

112  FORMAT!'  JFNUGTNX  ' ) 

5 0 1 R E A 0 ( b 1 T K A D 5 

» R I T E ( 6 * 1 0 b ) ( TRANS ( J ) ,J=1  , JFNU  1 
IT  X = N X - J r N U 


1 


IF (NX.gT • 1 1G0  To  50 G 
REAU(S,10y)  jfnu 
IF(JFNu.NF.O)g0  to  505 

502  continue 

1 F ( 1 LA .GE . IL AY  ) GO  TO  fi 7 p 
Go  TO  077 
8 9 0 W r J r E ( 6 . H 9 1 ) 

091  FORMAT!'  V 2 P GREATER  THAN  V2  - EKROi.-  STOP'l 

878  continue 
rewind  b 
202  call  exit 


77 

Format 

IE12,5,F7.2I 

8 1 

F 0 K M A T 

( 7 1 1 0 • 3 > 

8 5 

Format 

( 6 F l f)  . 3 ) 

79 

F 0 R M A T 

('  HRLSSuRl  = ' iE  1 2.5  TEMPER  AT,,hE  .'.F7.21 

83 

F URN  A T 

(3X, ' WATER'  iOXt,C02’» ft X,'OZONL'|7<,MT2o*»7X,'CO’*8A» 

1 * , »02» 

, HX  ) 

8 7 

format 

(•  Vl  ■ • ,F  10. 3,  • VZ  « ' .F  1 0. 3 , 'Ov»  ' ,F  1 1).  3 . • HOUND  «',Fl 

91 


■ 


‘ - ^ 


& <»*V- f A+i  r< 


t = • , r l o . 

3 , • (J  ( i.  V s'  ,1  1(j.  3 ) 

0 V £ o K i i A | 

( • p At'  I T Y L KF  Oil  E 

N c norm  ru  p 

A T • , F 17.31 

VI  * CROAT 

(•  Fill;  Of  flu.  fc  NCO'JNT  fc  A tO  • , 

IS) 

V 7 F 0 K M A r 

( * v l.l  0 T s’  |F  1 2 • 3 » 

• V I OP  = • ,F 1 2 . 3 , ’ SNU  = • 

,F 1 2. 3 , • I 

• . I ti  ) 

i d 3 format 

l • JMUJ  = • , i 1 

10  4 1 OKMAf 

lb<*  FlUoUt-MCY 

T P A N S • • 

) ) 

V 3 F o K M a r 

( ! 1 u ) 

VS  f o R M A T 

( f 1 U . 3 , E I 0 . 3 , F S . 3 

,F 1 0. 3 , 3bX  , 

13.14,13) 

loi  format 

l 1 b , 2£  1 □.  4 ) 

i o s format 

( b ( F 1 0 . 3 , E 1 2 . b ) ) 

Stop 

' t NO 

I3XQT 

3 7 

l .n 

I . 0 l . 0 

1 . 0 

t > o 

1 .0 

1 • 0 U 1 

l .uovs 

1 . U 1 9 2 1 . 0 3 ? 7 

1 .0602 

1 • o 7 1 9 

1 .0931 

1 . 12  6 9 

] .00  4 

I . I.IU7  1 . U l 3 

1.022 

1.033 

1.046 

1*066 

1.017 

1.03  1 . U ‘1  A 

1.072 

1 • 1 

1.127 

1.1/ 

1 .0 

1.0  t.o 

1 .0 

1 . o 

1 . G 

1 .0 

1 .11 

l.l)  1 .001 

l .00  2 

1 .Onv 

1.007 

1 . 0 1 1 

1 .0 

1.(1  1 .0 

1 . 13 

1 . o 

1 .0 

1 .001 

7*1, 

0 0 U (.  - 2 S 

6 S U • 

6 70  . 

.00  1 

1 b . 

. 

00b 

.U1 

1 6 0 

E -0  2 1 6 0 . 6 b 

. 100 

E-0  1 

180.6b 

o E 

0 6 . 7 8St  + 1 8 

.0 

£ 0 

. 0 £ 

0 

• n 

E 

0 

. 0 

E 

0 

.0 

E 

1 00 

E - 0 1 180.6b 

. 4 u 

E-Ol 

2 l 0. 22 

0 F 

0 6 . / a S L ♦ l 6 

• 0 

IT.  0 

. 0 £ 

0 

• i) 

F 

0 

.0 

E 

0 

.0 

E 

4 00 

E - 0 1 2 1 0 . 2 2 

.70 

e - o i 

22b .97 

0 E 

0 6«/Hbt+l8 

• c 

L o 

. 0 E 

0 

. n 

F. 

0 

. 0 

F. 

0 

.0 

E 

700 

E-Ot  22S.V7 

.10 

F Oo 

2 3 b . 6 S 

0 £ 

0 6 . 7 b b L * 1 8 

• 0 

E 0 

. 0 £ 

0 

* (■ 

£ 

0 

, 0 

E 

0 

.0 

E 

1 00 

E GO  2 3 b . 6 S 

.40 

£ Ou 

7 6 4 . lb 

(1  £ 

G 6.78SL-18 

. u 

11  o 

. a e 

0 

• n 

E 

0 

. 0 

E 

0 

.0 

E 

400 

E GO  2 6 4 . 1 S 

.70 

E 00 

270,6b 

0 t 

0 6.78SL+18 

. 0 

£ 0 

.0  E 

0 

• 6 

E 

0 

.0 

£ 

0 

.0 

E 

700 

E GO  270.6b 

. 10 

F ♦ 0 l 

270.6b 

0 E 

0 6 . 7 6 b t + 1 8 

• G 

£ n 

. 0 E 

0 

• 0 

E 

0 

.0 

E 

0 

.0 

E 

1 00 

E ♦ U l 2 7 0 . 6 S 

. 1 3 

E + 0 1 

267.17 

0 E 

0 6 . 7 8 b t + 1 6 

.0 

L 0 

.0  E 

0 

• 0 

E 

0 

.0 

F. 

0 

.0 

E 

1 30 

E ♦ 0 1 2 6 7 .17 

. 1 6 

£♦0  ) 

2 6 2 . S 1 

. 

0 E 

0 6 . 7 8 3 E ♦ 1 ft 

. (i 

£ 0 

. 0 F. 

0 

• c 

£ 

0 

. 0 

E 

u 

.0 

E 

160 

E ♦ 0 1 2 6 2 • b 1 

.20 

E + 0 l 

237.87 

0 t 

0 6 . 7 U S E ♦ 1 8 

.0 

E 0 

. U L 

0 

• r. 

F 

0 

.0 

E 

0 

.0 

t 

20 

£♦01  2 5 7 . t)  7 

, 2 S 

t ♦ 0 1 

2S3.20 

0 E 

0 6 . 7bbE+  l fl 

.0 

E 0 

.0  E 

0 

. n 

E 

0 

. 0 

E 

0 

.0 

£ 

2b 

E ♦ 0 1 2 S 3 . 2 0 

. 3 0 

£ ♦ 0 1 

249.36 

0 E 

0 6.78SE+I8 

• G 

L 0 

. 0 C 

n 

.0 

E 

0 

. 0 

E 

0 

.0 

E 

3 0 

£ + 01  24V. 38 

. 4 0 

£♦0  1 

243. 6b 

0 £ 

o 6 . 7 a b £ ♦ i a 

.0 

E 0 

.0  E 

0 

* n 

£ 

0 

.0 

t 

0 

.0 

E 

40 

£♦01  243.6b 

.so 

E + 0 1 

2 39 . 27 

0 E 

0 6.78SE+18 

. 0 

t.  o 

.0  £ 

0 

• 0 

e: 

0 

. 0 

E 

0 

.0 

£ 

SO 

£ ♦ U 1 239.27 

• 6 S 

c.  + 0 1 

2 3 4 , 0 9 

(J  t 

0 6 . 7 a S L ♦ 1 8 

. 0 

E 0 

. 0 E 

0 

. n 

e 

0 

. 0 

E 

0 

.0 

£ 

6b 

£♦01  2 3 M . u 9 

. ao 

£+0  J 

230.06 

0 K 

0 6. 7 c* St  ♦ 1 tt 

. 0 

I o 

. 0 £ 

0 

. n 

E 

0 

. 0 

F 

0 

.0 

E 

HO 

E ♦ 0 1 2 311.06 

. IU 

l.  *02 

227.06 

0 

o 

0 

u 

0 

0 

0 

u 

0 

0 

G 

G 

0 

0 

0 

u 
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I 


.0 

F. 

0 6 , 7 6 b t;  ♦ 1 8 

. n 

t.  0 

. o t: 

0 

. 0 

F 

0 

.0 

F, 

0 

.u 

E 

0 

. 1 u 

{*02  227.06 

. 1 3 

i +0? 

723.93 

. n 

E 

II  6 . 7 6',,  * I 8 

. 0 

1.  0 

. U F 

0 

• 1) 

f: 

0 

. 0 

t 

0 

.0 

L 

0 

. 1 J 

{.*1/2  2 2S  . VS 

. 1 6 

r ♦ 1/7 

;>  7 9 . 3 9 

. Ij 

E 

t)  6 . / llbl  ♦ 1 8 

. 0 

1 u 

• 0 1. 

0 

.0 

r 

0 

.0 

F 

u 

.0 

L 

0 

. 1 6 

t ♦ l*  2 2 28. 1.9 

. 20 

1.  + 0 2 

223.1? 

. u 

t 

0 6 i 7 6 b E * 1 8 

. U 

E 0 

• 0 1. 

0 
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